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1 Introduction and main results 



An example of a planar polygonal Markov field, referred to as the Arak process throughout 
this paper, was first introduced by Arak [1]. Further developments are due to Arak & 
Surgailis [2], [3], Surgaihs [16], Arak, Clifford & Surgaihs [4]. In this paper we focus our 
attention on polygonal Markov fields with V-shaped nodes, which are a particular class of 
ensembles of self-avoiding polygonal loops (contours) in the plane, interacting only by the 
requirement of disjointness. Not unexpectedly, these objects share a number of properties 
of the two-dimensional Ising model, including the presence of spontaneous magnetisation 
and absence of infinite contour nesting in low temperature region, see NichoUs [13] and 
Schreiber [15]. An important property of the Arak process and its length-interacting 
Gibbsian modifications is their isometry invariance. One might be tempted to view these 
purely continuum polygonal fields as a kind of continuum version of the Ising model. For 
low temperatures this opinion seems to be well founded. There is a number of relevant 
differences though in the much less understood high temperature region. In sharp contrast 
to the Ising model it is not clear how to define the infinite temperature non-interacting field, 
since some spatial correlation is always present due to the imposed polygonal nature of 
the contours. Therefore as the reference field for length-interacting Gibbsian modifications 
we choose the original Arak process, enjoying a number of striking properties including 
consistency, exact solubility and admitting the so-called dynamic representation in terms of 
equilibrium evolution of one-dimensional particle systems tracing the polygonal boundaries 
of the process in two-dimensional space-time, see Arak & Surgailis [2] and the Appendix 
below for details. 

The purpose of this paper is to show that, in analogy with the Ising model, the phase 
separation phenomenon is present for length-interacting polygonal Markov fields and it is 
gouverned by the Wulff construction, see Bodineau, loffe & Velenik [6] for an extensive 
reference. We establish our main Theorem 1.2 in the DKS set-up, as introduced by Do- 
brushin, Kotecky & Shlosman in their seminal monograph [7], and we only work at low 
enough temperatures. As a crucial tool replacing cluster expansion techniques and very 
well suited for our geometric setting we use a graphical construction built on contour birth 
and death process, as introduced by Fernandez, Ferrari & Garcia [9], [10], [11], see Sub- 
section 1.2. We took advantage of the particular properties of polygonal fields in order 
to characterise the model-specific surface tension, defined in Subsection 1.3, in terms of 



2 



hitting probabilities of appropriate planar random walks in random environment provided 
by the graphical construction. Even though we were only able to establish relatively weak 
results for the quality of approximation of the surface tension by its finite volume versions, 
we used the isometry invariance of the model to circumvent this problem. A particular 
feature of our approach is that rather than imposing periodic [as in Dobrushin, Kotecky 
& Shlosman [7]] or fixed sign boundary conditions [as loffe & Schonmann [12]], we work 
directly in the thermodynamic limit conditioned on the event that no large contours hit 
the boundary of the region. This allows us to avoid technical difficulties which would arise 
if we had to control our surface tension estimates in close vicinity of the boundary. Finally, 
the micro-canonical constraint considered in this paper requires that the excess of total 
magnetisation be larger or equal rather than just equal to a given positive threshold value 
- this avoids a number of technical complications which would otherwise arise due to the 
continuum nature of our setting and allows us to work with weaker versions of moderate 
deviation estimates and to rely upon general local central limit (LCL) results available in 
the literature rather than establishing an LCL theorem in its full strength specialised for 
our model. 

In analogy to the original DKS approach, the crucial ingredients of the proof of our 
main Theorem 1.2 are 

• the coarse graining estimates of Section 6, based on skeleton techniques slightly 
modified and specialised for our particular setting. The graphical construction of 
Subsection 1.2 is used as a crucial tool replacing cluster expansion techniques, 

• moderate deviation estimates for cut-off ensembles, stated in Section 3 and estab- 
hshed by the classical method of studying the restricted phase modified by actions 
of finely-tuned external magnetic fields, see e.g. Section 2 of loffe & Schonmann [12]. 
The graphical construction of Subsection 1.2 admits an extension for these modified 
fields. 

We believe the techniques developed in the present paper should in principle be appli- 
cable to general continuum models exhibiting isometry invariance and admitting polymer 
representation. 

The remaining part of the introductory section is organised as follows. Below, in Subsec- 
tion 1.1 we give a formal construction of the polygonal Markov fields. The next Subsection 
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1.2 is devoted to the graphical construction. The surface tension specific for our model is 
defined in Subsection 1.3. Finally, our main results are formulated in Subsection 1.4. 

Throughout the paper we make extensive use of the 'O, fi, ©' notation, with 0{X) and 
n{X) standing respectively for quantities bounded in their absolute value above and below 
by a constant times X, and with Q{X) — 0{X) D Moreover, we use c, C, Ci, C2 etc. 

to denote generic constants which can change their values from one statement to another. 

1.1 Length- interacting polygonal Markov fields 

The formal construction of the basic Arak process with empty boundary conditions in a 
bounded open set D C goes as briefly discussed below [wc refer the reader to [2] and [4] 
for further details]. In the sequel we assume that the boundary dD is piecewise smooth. 
We define the family Fd of admissible polygonal configurations on D by taking all the 
planar graphs 7 in D such that 

(PI) 7naD = 0, 

(P2) all the vertices of 7 are of degree 2, 

(P3) the edges of 7 do not intersect, 

(P4) no two edges of 7 are co-linear. 

In other words, 7 consists of a finite number of disjoint polygons fully contained in D and 
possibly nested. Further, for a finite collection (l) — of straight lines intersecting 

D we denote by F£)(Z) the family of admissible configurations 7 with the additional prop- 
erties that 7 C IJ"^^ and 7 fl is a single interval of a strictly positive length for each 
— 1, ...,n, possibly with some isolated points added. Let Ad he the restriction to D 
of a homogeneous Poisson line process A with intensity measure given by the standard 
isometry- invariant Lebesgue measure /i on the space of straight lines in IR^. One possible 
construction of /i goes by identifying a straight line I with the pair (0, p) e [0, tt) x M, where 
(psin(0), pcos(0)) is the vector orthogonal to I and joining it to the origin, and then by 
endowing the parameter space [0, tt) x M with the usual Lebesgue measure. With the above 
notation, the basic polygonal Arak process Ad on D arises as the Gibbsian modification 
of the process induced on F^) by A^), with the Hamiltonian given by the double total edge 
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length, that is to say 

IEy'^,cr (\ ^nr-expf— 2 length(7)) 

for all G C r^) Borel measurable, say with respect to the usual Hausdorff distance topology, 
and where r£)(A£)) denotes ^^{l) as defined above with I set to be the collection of all 
straight lines of A^- Note that by the total edge length length (7) of a polygonal configuration 
7 we mean here and below the sum of lengths of all constituent polygons. The expectations 
in (1.1) are taken with respect to the randomness of Ad- It should be mentioned at this 
point that in the literature on consistent polygonal fields one usually considers free rather 
than empty boundary conditions, see [2] and the Appendix below, yet the empty boundary 
object is better suited for the graphical construction below and for our further purposes. 

For a positive inverse temperature /3 > we consider the length-interacting Arak pro- 
cess ^ determined in distribution by 

exp(-/31ength(7)) 
dCiAo)^^^' Eexp (-/? length (^lij))' ^'^ 

with £(•) standing for the law of the argument random object. As shown in Theorem 3 
and Corollary 4 of [15] for (3 large enough the polygonal fields A^^i D C M."^, admit a 
unique whole plane thermodynamic limit without infinite contours, denoted in the sequel 
by A^\ see also below for its construction. The field A^^^ is isometry invariant. The 
thermodynamic limit A^^^ can also be shown to exist for /3 = 0, in the sequel it is denoted 
by A and its construction is given in the Appendix. 

It is known that for the inverse temperature /3 sufficiently large (in particular, for all 
(3 within the validity region of the graphical construction below) the thermodynamic limit 
j[[f3] exhibits only finite contour nesting, see Nicholls [13] and the discussion following 
Corollary 4 in Schreiber [15]. Whence, the contour ensemble ^[^1 partitions the plane 
into a unique infinite connected component {the ocean) and a countable number of finitely 
nested bounded regions (islands). We colour black and white the polygonal regions of 
this partition by declaring the infinite ocean white and by requiring that the collection of 
interfaces between black and white regions coincide with the collection of contours A^^\ 
which uniquely determines the colouring. Denote the resulting union of black regions by 
black [^[/^l] and the union of white regions by white [^[^1]. For a bounded region L'^ C let 
M[/ (^^^1) be the magnetisation in U determined by the coloured contour ensemble A^^ 
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under the assignment black i— > +, white i— > — . In other words, Mu [A^^^) is the total area 
of the black-coloured regions in U minus the total area of white-coloured regions in U : 

Mu {A^^) := Area (black[>l[^]] nU) - Area (white[>l[^l] n U) . 

For L > we shall abbreviate Mb.,(l) (A^^^) to Ml[P], where M2{L) stands for the disk of 
radius L centred at 0. The isometry invariance of the infinite- volume field A^^'^ implies that 

EMt; {A^f^^) = Area(C/) M[p], M[p] e (-1, 0), (1.3) 

where | M[(3] \ is further referred to as the specific spontaneous magnetisation at inverse 
temperature p. The 'black[-]' and 'M.(-)' notation will be also used for A^^^ replaced by 
a number of other polygonal fields enjoying the property that the corresponding contour 
ensemble determines a unique unbounded region, to be coloured white. 

1.2 Graphical construction 
1.2.1 Basic graphical construction 

As argued in Schreiber [15], the polygonal field admits a natural representation in 
terms of a graphical construction in the spirit of Fernandez, Ferrari & Garcia [9], [10], [11], 
which will be a crucial tool in our argument in the sequel, as replacing cluster expansion 
techniques. Below, we provide a description of this construction borrowed from [15]. Con- 
sider the space Cd consisting of all closed polygonal contours in D which do not touch the 
boundary dD. For a given finite configuration {I) :— {li, Z„) of straight lines intersecting 
D denote by Cd{1) the family of those polygonal contours in Cd which belong to Toil)- We 
define the so-called free contour measure ©£> on Cd by putting for C C. Cd measurable, 
say with respect to the Borel cr-field generated by the Hausdorff distance topology, 

Qd{C)= [ Yl exp(-21ength(^))(i/.*((/)) (1.4) 

with Fin(L[D]) standing for the for the family of finite line configurations intersecting D 
and where /i* is the measure on Fin(L[D]) given by dii*{{li, ...,ln)) '■— Y[i=i^l^{h) with /i 
defined in the discussion preceding (1.1). 

For /3 > we consider the exponential modification 0^' of the free measure Qd, 

e^\de) := exp(-/? length (^)) 9^3 (d^). (1.5) 
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Observe that for all bounded open sets D with piecewise smooth boundary the free contour 
measures ©£) as defined in (1.4) arise as the respective restrictions to Cd of the same 
measure © on C := U^i^(-n,n)2, in the sequel referred to as the infinite volume free 
contour measure. Indeed, this follows easily by the observation that restricted to 
coincides with for D2 C D^. In the same way we construct the infinite- volume 
Gibbs-modified measures ©'^1. The following result, which is Lemma 1 of [15] (note that 
the first result in this spirit is due to NichoUs [13], see Lemma in the Appendix ibidem) 
will be crucial for our further purposes. 

Lemma 1.1 For (5 >2 we have 

Q^\{e I dx n Vertices(^) ^ 0, length(^) > R}) < 87rexp(-[/3 - 2]R)dx, 

where the event {dxH Vertices(^^) 7^ 0} is to be understood that a vertex of 9 falls into dx. 
Moreover, there exists a constant e > such that, for /3 > 2, 

©[^1({^ I e Int^, length(^) > R}) < exp(-[/3 - 2 + e]R + o(R)), 

with Int 9 standing for the region enclosed by 9 ( recall that 9 E C is always a single bounded 
contour). 

Let V^ii3] be the Poisson point process on Cd with intensity measure ©^'. It follows then 
directly by (1.4) and by (1.1) that ^4^^ coincides in distribution with the union of contours 
in 'PqW conditioned on the event that they are disjoint so that 



( 

u ^ 



:i.6) 



where the conditioning makes sense because ^^{Cb) is finite as shown in Subsection 2.2 
of [15]. In particular, as argued in Subsection 2.2 and Theorem 2 ibidem, the law of 
is invariant and reversible with respect to the following contour birth and death dynamics 
(7f ) on r^. 

(C : birth[/?]) With intensity Q^{d9)ds do 

• Choose a new contour 6*, 
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• If 61 n 7f = 0, accept 9 and set 7^^^ := 7^ U 9, 

• Otherwise reject 9 and keep j^^^ig := 7^ , 

(C : death[/3]) With intensity 1 • ds for each contour 9 E remove 9 from 7^ setting 
7f 7f \ e- 

Moreover, £(^4^') is the unique invariant distribution of the above dynamics, see Theorem 
2 in [15]. These observations place us within the framework of the general contour birth 
and death graphical construction as developed by Fernandez, Ferrari & Garcia [9], [10], [11] 
and as briefly sketched below, see ibidem and Schreiber [15] for further details. Choose f3 
large enough, to be specified below. Define J^{C) to be the space of countable and locally 
finite collections of contours from C, with the local finiteness requirement meaning that 
at most a finite number of contours can hit a bounded subset of IR^. On the time-space 
R X J^{C) we construct the stationary free contour birth and death process {gs)seM. with 
the birth intensity measure given by O'^l and with the death intensity 1. Note that free 
means here that every new-born contour is accepted regardless of whether it hits the union 
of already existing contours or not, moreover we admit negative time here, letting s range 
through M rather than just M+. Observe also that we need the birth measure ©1^1 to be 
finite on the sets {6* e C | 6* n ^ 7^ 0} for all bounded Borel C in order to have the 
process (^s)seM well defined on R x J^{C). By Lemma 1.1 this is ensured whenever /3 > 2. 
To proceed, for the free process {gs)seM. we perform the following trimming procedure. 
We place a directed connection from each time-space instance of a contour showing up 
in (f?s)seR and denoted by ^ x [sq, Si), with 9 standing for the contour and [sq, Si) for its 
lifespan, to all time-space contour instances 9' x [s'q, s[) with 9'r\9 ^ (I), Sg < Sq and s'^ > sq. 
In other words, we connect 9 x [sq, Si) to those contour instances which may have affected 
the acceptance status of x [sq, Si) in the constrained contour birth and death dynamics 
(C) as discussed above. These directed connections give rise to directed ancestor chains 
of time-space contour instances, following [11] the union of all ancestor chains stemming 
from a given contour instance 9* = 9 x [sq, Si), including the instance itself, is referred to 
as its clan of ancestors and is denoted by An(^*). More generally, for a bounded region 
U in the plane we write Ans{U) for the union of ancestor clans of all contour instances 
9 X [sqi Si) with n f/ ^ and s G [so,Si). Lemma 1.1 allows us to apply the technique 
of domination by sub-critical branching processes, developed in [9], [10], [11], in order to 
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conclude that there exists Pg such that for each (5 > (5g there exists c := c(/3) > such that 

P(diamAn,(B2(x, 1)) > R) < exp(-ci?), s e R, x e R^ (1.7) 

with 182(3;, 1) standing for the radius 1 ball in centred at x. In the sequel we shall 
always assume that (3 > j3g, that is to say that (3 is in the validity region of the graphical 
construction. We see that for j3 > f3g all the ancestor clans are a.s. finite and we can 
uniquely determine the acceptance status of all their members: contour instances with no 
ancestors are a.s. accepted, which automatically and uniquely determines the acceptance 
status of all the remaining members of the clan by recursive application of the inter-contour 
exclusion rule. In this case, discarding the unaccepted contour instances leaves us with a 
time-space representation of a stationary evolution (7s)se]R on J-'{C), which is easily checked 
to evolve according to the whole- plane version of the dynamics (C) above. In Section 4 and 
Theorem 4 of [15] we argue that for all s e R the polygonal field 7s coincides in distribution 
with the thermodynamic limit (see Section 3 ibidem) for A^^^ without infinite contours, 
which is unique (see Corollary 4 ibidem). It should be observed that for each s e R the 
free field Qs coincides in distribution with the Poisson contour process Vqw . Since almost 
surely we have % Q Qs, we get the stochastic domination of the contour ensemble .4'^] by 

We also consider finite- volume versions of the above graphical construction, replacing 
the infinite- volume birth intensity measure O'^l with its finite- volume counterparts 0^' for 
bounded and open D with piecewise smooth boundary. Clearly, the graphical construction 
yields then a version of the finite- volume contour birth and death evolution (C). For each 
D denote by (7^ )se]R the resulting finite-volume stationary process on the space T{Cd) of 
finite contour configurations in D and write (^f )sg]R for the corresponding free process. It 
follows by Theorem 2 in [15] that 7^ coincides in distribution with for each s G R. 
Likewise, coincides in distribution with V^-m]. 

By representing the measures 0^^ as the corresponding restrictions of 0'^' we obtain 
a natural coupling of all the processes 7f,i?f,7s and Qs on a common probability space. 
We shall also consider coupled on the same probabihty space by putting = 7^. 
Likewise, we put A^^^ = 70, VqW] = Qo, 'T^nW = Qn ■ This coupling, referred to as the 
canonical coupling in the sequel, will be assumed without a further mention throughout 
this paper. 
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A simple yet useful application of this coupling is that 

Area(D) M[l3] - EMd{A^§) = 0(Area((9D © B2(l))) 

with © standing for the usual Minkowski addition [i.e. X ®Y :— {x + y \ x e X, ye Y}]. 
Indeed, this is immediately seen by observing that, by Lemma 1.1 and in view of (1.7) 
stating the exponential tail decay for ancestor clan diameters, under the canonical coupling 
of and A^^^ , the probability that a given point x e D is assigned different colours by 
these ensembles decays as exp(— Q(dist(x, Integrating over D and using (1.3) we 

obtain the required relation. 

1.2.2 Modifications of the basic graphical construction 

Below we discuss a number of modifications of the graphical construction, which will be 
of use for our further purposes. Apart from the area-interacting modifications all the 
remaining ones can be defined on the probability space of the basic construction, thus 
extending the canonical coupling. 

Imposing forbidden regions A particular property of the graphical construction which 
will be crucial for our further purposes is that it admits, on the same probability space, 
conditional versions on the events of the type no contour of the polygonal field hits [inter- 
sects] a given region U. Indeed, let U he a bounded subset of the plane R^. Then, adding 

the rule that all new-born contours hitting U [intersections of U only with the interior of a 
contour are not taken into account] be immediately discarded, to the trimming procedure 
constructing (7^) out of (qs) or, cquivalcntly, to the dynamics (C), wc obtain a stationary 
and reversible process {"fs-.u) easily seen to enjoy the property that the distribution of '-/g-.u 
for each fixed s coincides with the law of A^^^ conditioned on the event that A'^^^ DU = $. 
Put ^jgl^ := 7o:(7- Likewise, we define the conditioned version {"f^jj) of the finite volume 
process (7^ ) for which the distribution of coincides for each s G M with the law of 
.4.^' conditioned on {A^§ fl ?7 = 0}. We put A^j^ljj := 7,^^. In full analogy with the similar 
discussion above, the conditioned field -^jgi^^ is stochastically bounded by the Poisson con- 
tour process Vqw-u := {6 G V^m \ 6 nU = 0}. Likewise, A^^.jj is stochastically bounded 

Cut-ofF ensembles An important family of processes we embed into the original graph- 
ical construction are the cut-off ensembles for A^^^ They are defined as follows. For a 
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positive cut-off threshold a and a bounded region y C we consider the measure 
which is the restriction of O'^l to the family of polygonal contours which either do not hit 
V, or if they do hit V then their diameter does not exceed a. In this context, it is conve- 
nient to say that a contour 7 is a-large iff diam(7) > a and that it is a-small otherwise. 
Using bI/^I^"^'^ instead of O'^' for the contour birth intensity in the graphical construction 
we obtain a-cut-off version (7^ ^''*'^)sgR of the process (7s)seR (equivalently, we can simply 
reject all a-large contours hitting V upon their birth in the original graphical construction, 
which naturally extends the canonical coupling). Put ^t^l'"'^ := 7^^ It is easily seen 
that the a-cut-off polygonal field ^I/^l'"'^ := 7^^ '"'^ coincides in law with ^''^^ conditioned 
on the event that no contour hitting V has its diameter larger than a. Likewise, we con- 
sider with obvious definition the finite volume cut-off processes (7^ for open and 
bounded D with picccwise smooth boundary. Clearly, the finite volume a-cut-off polygo- 
nal field .4.^''"'^ := 7(f '"'^ arises as Ji^ conditioned on the event that no contour hitting 
V is a-large. In analogy with the similar discussion above, the cut-off field ^[^1'"'^ is 
stochastically dominated by VQ\i3]:a,v and ^^^'"'^ is stochastically dominated by P„[/3];c«,v. 

Clearly, we can combine the cut-off operation with imposing a forbidden region which 
leads to processes 7^;/"'^, 7^[/*'^, vA^i'^'^ and «4^!^'^ with obvious definitions, stochasti- 
cally dominated by Vq[s];cx,v.ij and Vf.[i3];a,v „ respectively. The canonical coupling is ex- 
tended in the obvious way. 

Area-interacting fields The final modification considered involves introducing an area- 
order term to the Hamiltonian of (1.2). To this end, for a bounded region 1^ C and 
/i e R we consider the polygonal field Ji^f^^ on W, given in distribution by 



Note that, unhke the field J}^ ^ is defined as a Gibbsian modification of the ther- 
modynamic limit ^[^1 restricted to W rather than as a Cibbsian modification of the finite 
volume field Aw- In particular, the laws of the fields Ai^f^ and Ai^ do not coincide; in 
fact Ai^^^ coincides in distribution with A^^"^ fl W. We will mainly use the area-interacting 
modification combined with the cut-off operation. The field ^[^''^^'"'^ is given in law by 




[7] 



cxp(/iM^y(7)) 



dC{A^P^ n W) 



Eexp(/iMiy(^['3]))' 




[7] 



exp(/iMH/(7)) 



(1.8) 



C^£(^[/3];a,y n W) 



Eexp(/iM^(^[/3]i«'^))' 
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To proceed with the graphical construction we assume that 




(1.9) 



and observe that adding a single a-small contour to a contour configuration 7, 7 fl ^ = 
0, can change the magnetisation Mvk(7) by at most tt length(^)^/2 whence the value of 
hM.\.Y{j) can change by at most (3\ength{6)/2. With 7 standing for the current contour 
configuration, we modify the original graphical construction by 

• constructing the free birth and death process Qs, s G R, with birth intensity measure 
Qimh^y and death intensity 1, 

• at the trimming stage, by accepting a time-space contour instance 9 x [sq, si) 

— with probability if hits 6' for some previously accepted contour instance 
9' X [sq, s'l) alive at time Sq, 

- with probability exp (-| length(6') + h[Mw{'y D 9) - Mwi'f)]) otherwise. 

Observe that the last probability falls into (0,1] because of (1-9). Denote the resulting 
trimmed process by 7^. The validity of this construction requires a justification. In fact, 
we have to redefine here the notion of an ancestor clan. We set a directed connection from 
a contour instance 9* = 9 x [so,Si) to all contour instances 9'* = 9' x [sq,s'i) such that 
Int ^^nlnt 9' ^ ^ (which is weaker than the condition 9r\9' 7^ of the original definition) and 
s'o < So I s'l > So- Clearly, these are all contour instances which may affect the acceptance 
status of 9*. The union of all the directed chains stemming from 9* is called the ancestor 
clan of 9* and denoted by An(^^*). Likewise, for s G M and C we write Ans{U) for 
the union of all the ancestor chains of contour instances 9 x [sq, Si) alive at time s (i.e. 
So < s < Si) and such that Int^flC/ 7^ 0. In full analogy with (1.7), Lemma 1.1 guarantees 
that for P large enough (larger than some Pg) we have 



for some c = c{(3) > 0. Clearly, this imphes that the ancestor clans are a.s. finite, thus 



so that (1.10) holds. It follows by the general theory developed by Fernandez, Ferrari & 
Garcia [9], [10], [11], that so constructed % for each fixed s coincides in law with 



P(diamAn,(B2(a;, 1)) > R) < exp{-cR), s eR, xeR^, 



(1.10) 



ensuring the validity of the construction. In the sequel we shall always assume that f3 > $g 
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Moreover, it is easily seen that, for each s e R, ^5 coincides in law with the Poisson contour 
process VQii3/2];ot,vr\W. Consequently, the almost sure inclusion % C yields the stochastic 
domination of A^^^'""'^ by Vqi0/2];c,,v fl W. 

Clearly, the above construction can be easily extended to take into account forbidden 
regions. For bounded measurable [/ C denote by the polygonal fields arising 

by conditioning ^[^''^''"'^ on none of its contours hitting U. It is easily seen that 
can be represented by the graphical construction of this paragraph, with the additional rule 
that all contours hitting U be immediately discarded. In analogy with a similar observation 
made above for ^{^'"'^, also here it should be noted that A^^'^ij"'^ coincides in law with 
A^l:^'^ n W rather than with A^^^'^ . 

Moreover, in full analogy with the argument above, wc sec that A^^'.^'^ is stochastically 
dominated in the sense of inclusion by the Poisson contour process 'P0[/3/2];a,v.^ fl W. 

Note that the above construction provides a natural coupling for area-interacting fields 
with cut-off and (possibly) forbidden regions imposed, under the constraint (1.9). To 
distinguish it from the canonical coupling available for fields with no area interaction as 
discussed above, we shall call this coupling the canonical coupling for area-interacting fields. 

1.3 Surface tension 

The purpose of this section is to define the surface tension functional specific for our model. 
To this end, for a given bounded and convex domain D CM.^ and S > we consider the 
family C^^'^ of self-avoiding polygonal paths in D connecting the balls B2(a;, S) C. D and 
182(1/, 5) ^ with the additional property that the first and last segments of the paths 
do not intersect the interiors of the balls M2{x,S) and M2{y,S) respectively, but they do 
touch their respective boundaries and the intersection points coincide with the initial and 
final point of the path. In other words, moving along a path in C^^'^ we travel from 
dM2{x,6) to dM2{y,S), with the initial segment falling outside 182 (x, 5) and with the final 
segment outside 182 (y, S), which does not prevent us though from passing through 182(0:, S) 
and 182(1/, ^) along the remaining segments. Next, we introduce on C^^'^ the free measure 
Q'd^^'^, constructed in full analogy with the definition of the free contour measure as given 
in (1.4). For a finite configuration (/) of straight lines crossing D write C^^'^{1) for the 
collection of those paths in C^^'^ which only contain segments of the lines in (/) and 
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exactly one non-zero length segment on each line. For measurable C C C^^'^ we put 

er''\C)^ [ Yl exp(-21ength(e))d/.*((0) (1.11) 

with Fm(L[D]) and fx* as in (1.4). Note that the initial point and the endpoint of the 
path 9 in the above definition are uniquely determined, respectively as the intersection of 
the first and last segment of the path with dM2{x, S) and dM2{y, S). Likewise, we define the 
/5-tilted measures [0^*^^'^]['^1 by 

[Q^^y,Sf](^d0) :=exp(-/31ength(^))e^"^^''^(de). (1.12) 

As observed above for the free contour measures, also the path measures are consistent in 
that Q^^'^ — [Q^j^^'%D for D C D' and, consequently, we can construct the whole plane 
free measure Q^^V'^ and its tilted version [e^^y'^jW^ both defined on C^^^'^'^ = C^^^'^ := 

Uoo /^a;•f-»J/;<5 
n=l'^(-n,n)2- 

To proceed, write 

tIs)[x ^ y] / ^A^^ ne = iD)[Q''^y'Y\de). (1.13) 

Put 63; := (1, 0), fix some small 5 > and let 

rf :=-ilogT;[^l[O^Ae.]. (1.14) 
The surface tension is defined as the hmit 

rl^l := lim rf = - lim ilogr/f.^[0 ^ Ae^]. (1.15) 

A— »oo A— »oo A ^ ' 

It is clear that the asymptotic behaviour of rf as A — ^ oo does not depend on the choice 
of 5 above - indeed, changing 5 is easily seen to result only in a bounded and uniformly 
non-zero prefactor before T'^'f)[0 ^ Ae^.], which is negligible in the logarithmic large A 
asymptotics. This is why our notation does not take into account the dependency of rf^ 
on 5. The existence, finiteness and strict positivity of the limit in (1.15) for large enough 
and other properties of the surface tension are discussed in Section 5 below, see Lemma 
5.3 there. 
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1.4 Main results 



Our main result below states that, at low enough temperatures, conditioning the white- 
dominated phase of polygonal Markov field to contain black- coloured regions of total area 
exceeding its expectation by an area-order quantity results in aggregation of the excess 
black area and in formation of a macroscopic-size disk-shaped region (Wulff crystal) of 
black-dominated phase, separated from the outside white phase by a single large contour. 
Moreover, the probability of such area-order black exceedances exhibits perimeter-order 
exponential decay. 

As shown in Section 2 below, for a ^ log L with overwhelming probability there are no 
CK-large contours of A^^^ in B2(-L). Thus, it is natural to consider the regions separated by 
Jl(log L)-large contours of ^'^1 and to assign them, in the obvious way, black or white phase 
labels. In this language, we show in this paper that the single large contour determining 
the Wulff shape encloses a disk-shaped portion of black-labeled phase region surrounded 
by ocean of white-labeled phase. 

As already discussed in the introductory section, since our main results are formulated 
directly under boundary conditions induced by the thermodynamic limit ^'^1 rather than 
with periodic or monochromatic boundary conditions, we have to explicitly rule out the sit- 
uation where the phase separating curves cross or go along the boundary of the considered 
finite volume region B2(L). To this end, for a > shall write A/"[q;, L] for the event that no 
a-large contour of A}^^ gets closer than at the distance Go; to the circle Si(L) := dW)2{L). 
In what follows we shall write 



a 



[L] := v^logL. 



(1.16) 



Our main result is the following theorem. 



Theorem 1.2 For < a < 27r| M[/3]| we have 



P (Ml {A^I^^) > M[/?]7rL2 + aL\ J\f[a[L],L] holds) = exp 




Lr% + 0{a[L])\ 
(1-W) 




Moreover, there exists a constant Cia^ 



such that on the event 



{Ml (^I'^l) > M[(3]7rL^ + aL^, J^[a[L], L] holds} , 
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for sufficiently large L we have with probability arbitrarily close to 1 

• There is exactly one Ciarge log L-Zar^'e contour ^laxge, 

• This 6'iarge Satisfies 

with ph standing for the usual Hausdorff distance. 

Note that in the sequel we shall refer to the condition Ml (^A^^^) > M[/3]7rL^ + aL^ as to 
the micro- canonical constraint. 

The remaining part of the paper is the proof of Theorem 1.2 and is organised as follows. 
In Section 2 below we establish upper bounds on occurrence probabilities of large contours. 
Next, in Section 3 we study moderate deviation probabilities for cut-off contour ensembles 
of polygonal fields. Section 4 provides a simple yet important lemma allowing us to factorise 
the avoidance probabilities of A^^^ over disjoint regions. This is followed by Section 5 deahng 
with properties of the surface tension, and then by Section 6, where we establish coarse- 
graining skeleton estimates. The complementary lower bounds for occurrence probabilities 
of large contours are stated in Section 7. Finally, in Section 8 we complete the proof of 
Theorem 1.2 by putting together the results of previous sections. 

2 Exponential tightness bounds 

In this section we show that although the total length of the contour ensemble A^^^riB2{L) is 
clearly of the area order ©(-L^), this is due to the contributions of 0(logL)-small contours, 
while the contribution of Q(logL)-large contours is of order 0(1) with the corresponding 
large deviation probabilities exhibiting exponential decay. To put it in formal terms, with 
a > not necessarily given by (1.16), write Lq,l := Lq,l [A^^^] for the family of a-large 
contours of A^^ hitting M2{L) and, in general, let ]Lq,_l[7] stand for the family of ct-large 
contours of a contour collection 7 which hit B2(-L). We claim that 

Lemma 2.1 For each k < P/2 — 2 there exist M^C — C{P,k) < 00 such that for all 
a > ClogL and X> we have 

P (length(Lc,,L) > A) < Mexp(-«;A) 
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and the same applies for 'La,L replaced with Lq_l[A'^''"' ], La,L[vA'(^).'")'' ■*] and'ha,L[^.)^^'°''^ \ 
'^a,L[-^^.y^!^'^'\ for h within the validity range of (1.9). 

Note that it is natural to regard this lemma as an exponential tightness statement for 
length (Lq,^l), whence the title of the section. 

Proof We provide the proof for the polygonal field J^^'^ only, since the argument goes 
exactly along the same lines for all the modified fields obtained from the variants of the 
graphical construction discussed in Subsubsection 1.2.2 and admitting stochastically dom- 
inating Poisson contour processes. Note that the assumption n<(3/2 — 2 was imposed for 
the purpose of dealing with area-interacting processes with the external field h within the 
validity range of (1.9), which admit stochastic domination by the Poisson contour process 
Vq[i3/2] . For the remaining polygonal fields considered in the statement of the lemma, with 
no area interaction, a stronger stochastic domination by Vqw is available and the assertion 
of the lemma still holds if we choose k < (3 — 2 rather than k < (3/2 — 2. 

To proceed, use the graphical construction to conclude that the total length of contours 
in Lq,,l is stochastically bounded by the total length of a-large contours of Vq{i3/2] hitting 
l^a,L. Thus, by the definition of a Poisson point process. 



Eexp («;length(LQ,^L)) < exp 



{eec I enB2(L)5^0, e is a-large } 



< 



/■oo 

e''"C(«) + / e''\(A)(iA, 

J a 

where the last inequality follows by integration by parts with 

((A) := Q^'^'^\{e e C I enn^iL) ^ 0, length(^) > A}) 
In view of Lemma 1.1 this means that 

Eexp («;length(Lo,^i)) < exp 

with some constant Ci(/5, k). Thus, using Markov inequality we get 

Ci{P,K)L'^&cp{[K + 2-^]a) 



C,{(3,K)L^exp{[K + 2-(^]a) 



^/•, , /-.r V ,v Eexp (k IcngthfLo, r)) 

P length L,,i > A < ^ , ' < exp 

exp(/tA) 



which completes the proof of the lemma for a > ClogL with large enough C — C{(3,k). 
□ 
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3 Moderate deviations for cut-off ensembles 



The current section deals with the properties of the cut-off ensembles A^^'"''^ arising by 
conditioning the original field A^^^ on containing no cc-large contours hitting F C R^. 
Recall that we assume here that /3 > /3g and P > Pg so that P falls into the validity regions 
of the graphical construction discussed in Subsection 1.2 as well as of its area- interacting 
modification discussed in Subsubsection 1.2.2. We consider a not necessarily given by 
(1.16). Our first observation is that imposing a cut-off does not change the expected 
magnetisation by too much 

\7rL^M[p] -EMl I = 0(L^exp(-ca)) (3.18) 

with some c > 0. Indeed, this follows by the fact that, in view the stochastic domination of 
^[^1 by Vem and in view of Lemma 1.1, an a-large contour shows up in A^^^ fl B2(i>) with 
probability 0(L^ exp(— ca)), whence conditioning on the absence of this event can change 
the probabilities of other events by at most 0(L^ exp(— ca)), consequently the variational 

distance between the laws C^A^'^^) and £(^[^l'"'*2(i)^ jg of the same order 0{L^ exp(— ca)). 
To get (3.18) it suffices now to observe that the magnetisation over M2{L) is a.s. bounded 
in absolute value by nL^. 

Another useful observation is that the impact of imposing a forbidden region for cut-off 
ensembles can also be very well controlled. In formal terms, we claim that for a collection 
7 of a-large contours, a > 1, in B2(I/) we have 

EMl -EMl f^i^:!'"'^^^) = O (Area(7 e B2(l))) = O (length(7)) . 

^ ^ (3.19) 

This is an immediate consequence of the fact that, by (1.7), under the canonical coupling 

of g^Yid AI^I'"'^^^^'' the probability that the colours assigned to a given point x 

by these ensembles differ, is of order 0(exp(— cdist(a;, 7))), c > 0. 

The argument leading to (3.18) and (3.19) above can be easily modified to yield the 
following combination of these relations. Let 7 be a collection of a-large contours in M2{L). 
Then 

E Ml U 7) - I M[/3] I MM | = 0(length(7)) (3.20) 

provided a > C log L for sufficiently large C. 

The main result of this section is the following moderate deviation bound 
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Theorem 3.1 For each (3 large enough there exists a positive constant Ci = C'i(/3) such 
that, uniformly in L,a > Ci log L and in finite collections 7 of polygonal contours in 
we have for allO < A< Cf ^L^/logL 



P {A^t/^^'^ U 7) - EM^ ^^^^ U 7) I > a) < exp 

with some constant c > 0. 



a 



Proof Write 

/^L,.-M,(41-WU7) 

and let 

For /i e R consider the partition function 

Z[h] :=Eexp (/i/x^J 

The following estimate, vahd for all h satisfying (1.9), is the crucial ingredient of our proof: 

log Z[h] < /iE/x° ^ + /i^LVV2 (3.21) 

for some a > 0, uniformly in L,'y,a and h within the validity region of (1.9). To see that 
(3.21) suffices to complete the proof of the theorem, take first < A < ^^2^, put h :— 
which clearly satisfies (1.9), and then use Markov's inequality to conclude that 

P > E,i„ + A) < < eM>^'LV/2 - Ak) = exp {-^)^^^ 

Next, for A > choose k < /3/2 - 2 and C((3, k) as in Lemma 2.1 and assume that 

Ci in the present lemma is chosen so that a :— 2C(/3, k) logL < a. Then, on the event 
{a*l,7 > E//2,-y + A} there are two possible scenarios: 

• The total length of a-large contours in ha^L \ A!^}'.'^'^'^^^^ j exceeds which 



can 



happen with probability at most Mexp(— ^^) by the exponential tightness Lemma 
2.1, 
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• The total length of a-large contours in L^,!, ^^4^2!"''^^''^''^ does not exceed Since 
is the lower bound for the length-to-enclosed-area ratio for an a-small contour, this 
means in particular that the total area enclosed by a-large contours of A^l'"'^'^^^'' 
falls below A/ 4. Denoting by 7 the family of contom's ha^L (^-^j^j;"'*^*"^^ j U 7 and 
taking into account that the change of magnetisation induced by adding a contour is 
bounded in absolute value by twice the area it encloses, conditionally on given 7, we 
are led to 

E/.f < EMl J^'^-"^ U 7) + A/2. 
Now, in full analogy with (3.19), on the considered event we get 

EM, u 7) - = O (^) . 

Next, in full analogy with (3.18), we have 

l^f^l^- fit J ^0{L'eM-c&)) 

which goes to faster than the inverse of any polynomial under appropriate choice 
of C(/9, k) in Lemma 2.1. Putting the above relations together we conclude that 

Epil^<Epil^ + A/2{l + o{l)). 

Recalling that /i^^ given 7 coincides in law with /i2,7) we are led to 

P (/x2,, > E2,,|7) < P (/x!,^ > Ef,^ + A/2{1 - 0(1))) . 

Now, choosing Ci so that A < ^^f^, we can bound above the last probability by 
exp(— g^^a ) applying (3.22) with a and 7 replaced there by a and 7 respectively. 

Combining the above two points with (3.22) and noting that the probability of {fJ'l,^ < 
E/i2 ,y — A} can be dealt with in a completely analogous way, we obtain the assertion of 
the theorem. 

It remains to verify (3.21) for h satisfying (1.9). We extend the notation by putting 
Noting that 

^ log Z[h] = Enl'^;, ^ log Z[h] = Var[/.a (3-23) 
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and Taylor expanding the logarithm of the partition function up to the second order term 
yields 

log Z[h] = hE^,l^ + Y^^KoJ (324) 

for some h* between and h. Wc claim that, uniformly in L,a,'y and h satisfying (1.9), 
the variance Var[/x2'!|;] is of the area order 0{L'^) 

VarK;!;] = 0(L^), (3.25) 

which, once established, will immediately yield the required relation (3.21) as a conclusion 
of (3.24). To prove (3.25) we show that, for Ui,U2 Q M2{L), uniformly in 7 and in h 
satisfying (1.9) 

Gov [M^,(<(^^;?f U7);M^,(<fi-f U7)_ 
O (Area(?7i) Area(;72)[Area(t/i © B2(l)) + Area(t/2 © ©2(1))] e'^'^''^^^''^^^) 

(3.26) 

for a positive constant c, with © standing for the usual Minkowski addition. Indeed, with 
the representation provided by the graphical construction for area-interacting fields in Sub- 
subsection 1.2.2, conditionally on the event {Ano(f/i) C Ui © Ba C^''*^^^'^'^ ), Ano(f/2) C 

U, © B2(^ii^%^)} the random variables MuM^^y.^^'^'^^ ^ ^) ^uM^^^iy.^^'^'^^ U 7) 
are independent. But in view of (1.10) the probability of this event does not fall below 
1 - 0([Area(f/i ©12(1)) + Arca([/2 ©182(1))] exp(-cdist([/i,[/2))). This observation com- 
bined with the fact that \MuMB;'/^iy.:^''^^^)\ < Area(f/i), i = 1,2 implies (3.26). The 
required relation (3.25) follows now from (3.26) by usual argument based on splitting 
B2(-L) into ©(L^) disjoint regions of diameter and area 0(1) and then noting that, with 
the magnetisation contributions coming from distant regions exhibiting exponentially de- 
caying covariances, the asymptotic order of the total magnetisation variance Var[/i^'^] is 
determined by the sum of covariances between regions within distance 0(1) from each 
other, which yields the desired order 0{L'^). The proof is complete. □ 



Remark 3.2 We note that the bounds in Theorem 3.1 are of optimal order only for the 
probabilities of positive deviations {fJ^I^ > E,fi1 ^ + A} , A > 0. We believe that the probabil- 
ities of negative moderate deviations {a*2,7 < ■'^A*l,7 ~ ^} ^^^^ {i^l < E/xl — A}, A <^ 
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L^, exhibit Gaussian-type decay exjp{—D,{A^ / L^)) as in classical moderate deviation regime, 
in full analogy with similar phenomenon for the Ising model, see [8], (1.1.2), (2.3.2) in 
[12] or Section III.C.l in [6] and the references therein. Since this falls beyond the conteoct 
of our further argument, we do not discuss this issue in the present paper. 

As an easy corollary from Theorem 3.1 we conclude that 

Corollciry 3.3 With ^4 > 1, a and 7 as in Theorem 3.1 and with Ci in Theorem 3.1 large 
enough we have uniformly 



P (|Mi (^g^^ U 7) - EMi U 7) I > a) < exp (^-c 

with some constant c > 0. 



A' A 



VO(L^exp(-CQ;)) 



Proof This is a direct conclusion of Theorem 3.1 combined with the observation that 
the variational distance between the laws C{A^l.^ Pi M2{L)) and C{A^2.'^'^'^^^^ fl M2{L)) is 
of order 0(L^ exp(— cct)), in full analogy with the argument leading to (3.18) above. □ 

Another useful corollary relies on a straightforward observation that the proof of Theo- 
rem 3.1 apphes with only minor modifications for ^[/^l^^'^^C-t-) replaced by A^^'^j^^'^'^^^'^ with 
\h\ < H/a, H small enough. In formal terms, 

Corollary 3.4 With H > small enough, for each (3 large enough there exists a posi- 
tive constant C = C{(3,H) such that, uniformly in L,a > ClogL, finite collection 7 of 
polygonal contours in and \h\ < H/a, we have for all < A < C'^L"^/ logL 



'^jim^ML, u,)-EU, (.41";*-*' U -,) I > a) < exp 
with some constant c > 0. 



A^ A 
a 



We omit the proof of this corollary which is just a simple repetition of the proof of Theorem 
3.1. 

Below, we provide some further auxiliary results related to moderate deviation proba- 
bilities for cut-off ensembles. Note first that we can establish a bound analogous to (3.25) 
for the third cumulant of /x^'^ : 

^3 



log Z[h] = 0{L') (3.27) 
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uniformly in ct, L, finite contour collection 7 and h satisfying (1.9). We omit the details of 
the argument, based on the relation (1.10), since it goes along the same lines as the proof 
of Lemma 5.3 in Baryshnikov & Yukich [5] (valid for arbitrary order cumulants in fact). In 
particular, in view of (3.23) combined with (3.27), we get for h within the validity range 
of (1.9) 

E/.^;!; = Ei^l^ + /i Var[/.2,^] + 0(h'L') (3.28) 

uniformly in I/,7. To proceed, assume that 7 is a finite contour collection in M2{L) with 
Area(7 © B2(log^L)) < L^/logL. We claim that under this condition we have (3.25) 
strengthened to 

Var[/.2,g = e{L') (3.29) 

uniformly in 7, L. Indeed, observe that by (1.7) the probability of the event {Ano(7) ^ 
7 © B2(log^ L)} is of order at most 0(L^ exp(— clog^ L)), c > 0, whence under the canoni- 
cal coupfing with probability 1 — 0(L^ exp(— clog^ L)) the field .4^2.'"''^^*^^^ coincides with 
over the whole complement of 7 © B2(log^ L). Consequently, 

|Var[/i2]-Var[//2,^]| = o(L^). 

Now, mimicking the proof of (3.18) we check that 

I Var[/x2] - Var[/Xi]| = 0{L^ exp{-ca)) = o{L^) 

provided a > ClogL with C large enough. This will yield the required relation (3.29) as 
soon as we show that for the field A^^^ the variance of magnetisation has the required order 

Var[//L] =Q(L2). (3.30) 

To this end, we fix some large A > 0, large A; e N and small e > and we note that 

inf |var[Mz)(^gl)] | Area(L') > eA^ card Vertices (D) < k, \ength{dD) e [eX,e-^X]} > 0, 
^ \3.3i) 

with the infimum taken over all bounded domains D with polygonal boundary, possibly 
chopped off by intersecting with M2{L). Indeed, this can be proven by observing first that 
the mapping D 1— > (j){D) := Var[M£)(^^')] admits only strictly positive values and it is 
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continuous with respect to the pseudo-metric p'^{Di, D2) '■— inf^g]R2 ph{Di, x + 02). Thus, 
putting V[\k,e\ := {D \ Are&{D) > eA^, card Vertices(D) < k, length(9D) e [eA,e-U]} 
and noting that T>[X, k, e] is compact in p^, we see that 0o '■— infD6i5[A,fc,e] <P{D) > 0, which 
yields the required relation (3.31). To proceed, note that the variance Varf/x^] in (3.30) is 
bounded below by the expectation of the conditional variance of /il given the ensemble of 
external (outermost) contours 9 in A^^^ fl M2{L) satisfying the constraints of the infimum 
in (3.31) for D := IntO. Thus, taking into account that given the presence of such 9 the 
behaviour of the process A^'^^ inside 9 is independent of that outside 9 and then using (3.31) 
to conclude that each such 9 present brings a contribution of at least 0o to the considered 
conditional variance, we have Var[/i2,] bounded below by 0o times the expected number 
of external (outermost) contours 9 in A'^^^ fl M2{L) as in (3.31) with D = lnt9. Since this 
number is clearly of the area order Q{L^), the required relation (3.30) has been established, 
which completes the argument for (3.29). 

Putting together (3.28), (3.29) and the observations that h = o(l) by (1.9) and that 
E/x^'!|; strictly increases with h we come to 

Corollciry 3.5 There are positive constants Kq — Kq{[3) and C — C{(3) such that for 
each a > ClogL, each A with \A\ < K^L'^/a and each finite contour collection 7 with 
Area(7 ©B2(log^L)) < L^/logL there exists a unique value h — /t[A,L, 7] of external 
magnetic field such that 

and 

h = e(A/L2) 

uniformly in a, A, L, 7. 

Our next statement provides a lower bound for moderate deviation probabilities of A*2,7' 
complementary to the upper bound of Theorem 3.1. 

Lemma 3.6 For all < A <^ L'^/o:, with a and 7 as in Corollary 3.5 and with a <^ L/ 
logL we have uniformly in A, a, L, 7 

P {i^l^ > Ei^l^ + A) > cxp(-0([A + L\ogL]yL')). 
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Proof Write using Corollary 3.5, putting for brevity h[-] :— h[-, L, 7], 



PK,^ > Epil^ + A) > F{\pil^ - Epil^ - A - LlogLl < LlogL) > 
exp {-h[A + L log L] [Ei^l^ + A + 2L log L]) E exp {h[A + L log L]i^l^^) 

^naM^+Llo^L] _^^aMA+LlogL]\ ^ ^ 

and use Jensen's inequality to bound it below by 

p j^|^a,MA+LlogLl _ ^^a,MA+LlogLl| ^ ^ exp(-/i[A + L log L] [ A + 2LlogL]). 

Thus, taking into account that h[A + LlogL] = 6([A + L log L]/L^) by Corollary 3.5 
and that f^\^'^^M^+LiosL] _ ^^a,MA+LiogL]| ^ ^^^g^^ ^ ^ _ ^^^^ ^-^^ Corollary 3.4, 

completes the proof of the lemma. □ 

4 Decoupling lemma 

The purpose of this section is to establish Lemma 4.1 stating that the avoidance probabili- 
ties for the field over disjoint regions can be very well approximated by the product of 
the corresponding avoidance probabilities for individual regions. Even though this lemma 
is a direct conclusion from the graphical construction, we state it in a separate section due 
to its importance in our further argument. 

Lemma 4.1 Assume that Ui,U2, ■■■,Uk, k > 1 are disjoint bounded regions in such 
that miuj^j dist([/i, Uj) > A log[A; supf^^^ diam(t/j)]. Then, for some C > we have 

(k N 
A^f^^ n IJ = 

1 + O ( (log k) e-^^ ^ Area(t/i © Ball)) ) ) H ^ {^^^ n t/,- = 0) . 

Proof The exponential decay of ancestor clan sizes in the graphical construction (1.7) 
yields 

F{El) = 0(Area(C/i © B2(l)) exp(-CA)), i = 1, . . . , A; (4.32) 
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with 



8r-^ {Ano(C/,) CC/,eB2(A/2)}. 



Write Xj, j = 1, /c for the event 

X, {^[^1 n U, = 0} 

and use the canonical couphng of the graphical construction for J\^^ with the conditional 
graphical construction for the field .^u/oi as provided in Section 1.2 to conclude that 



\k/2\ 

i=Lfc/2j+l 




(4.33) 



-CA 



\k/2\ 

Lfc/2j+i /4.34) 



i=l 



Combining (4.32) with (4.33) leads to 

P^Qx.j = 1^1 + |^^Area(C/,©B2(l))e 
The assertion of the lemma follows now by recursive apphcation of (4.34). □ 

5 Existence and properties of surface tension 

This section deals with the existence of the hmit (1.15) defining the surface tension func- 
tional specific for our model. The argument below relies on a number of technical properties 
of the quantity T^^ [• <-> •] and is split into several subsections. Our main tool here is the 
random walk representation of surface tension, stated in Lemma 5.1, and our main effort 
is concentrated on establishing the crucial finite volume approximation Lemma 5.2. As 
everywhere in this paper, the results below are valid for j3 large enough. 



5.1 Optimising and freezing initial segments 

It will be convenient for our further purposes to switch between several alternative but 
asymptotically equivalent variants and representations of the surface tension. In this sub- 
section we argue that modifying and freezing the directions of the initial segments of the 
polygonal path in the original definition (1.13) of the functional T^^}^[x ^ y] does not alter 
its asymptotic behaviour for large dist{x,y). To this end we consider a version Tj^^[x ^ y] 
of the quantity Tj^^[x ^ y], which arises as the supremum of the integrals as in (1.13), 
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but in which the initial point of the first segment is now allowed in the whole M2{x, 5) 
rather than just on 9B2(x, 5), the endpoint of the last segment is allowed in the whole 
B2(y, S) rather than just on dM2{y, S), and the directions of both segments are fixed so that 
the integration is carried out over the remaining segments only. It is easily checked that, 
provided the distance between x and y is large enough, 

C-'tI^^ [x^y]< fl^j [x^y]< CT^^ [x ^ y] (5.35) 

for some C — C{P,S) > 1 independent of D,x,y. Indeed, the impact of taking the first 
and last segments fixed in the optimal way rather than integrating over them is easily seen 
to be only confined to close neighbourhoods of the initial point and the endpoint of the 
path, and can be compensated at a constant probability cost by appropriately adjusting a 
small number of initial and final segments. We also consider finite volume versions T^/^^j~^-^ 
and f(f^'^) of T(J^^ and f^], putting in analogy with (1.13) 

and defining T^^}^^^ in the same way as T^^}^ with the additional requirement that the 
whole path be contained in D. If the domain D contains neighbourhoods of x and y 
(say, B2(x, 25) C D and B2(y, 25) C D), a relation analogous to (5.35) is easily verified to 
hold for X and y far enough 

C-'Tj^U [^^y]< tI^U - ?/] < CTj^^D) - y] (5-36) 
with some C := C{j3,5) > 1 independent of D,x,y. 

We close this subsection with one more quantity, to be of use in the sequel, for which 
a relation analogous to (5.35) and (5.36) is valid. Write 

df [x ^ y] = / [e'^'^y-'Y^de) = [e^-2/;'5][^i(c^-2/;5) (5.37) 

and, as in the definition of f (';][• ^ •], let [x <-> y] be the supremum of integrals 
as in (5.37), but with the initial point of the first segment now allowed in the whole 
B2(x, 5) rather than just on dM2{x,5), the endpoint of the last segment allowed in the 
whole B2(y,5) rather than just on 9B2(y,5), and the directions of both segments fixed 
so that the integration is carried out over the remaining segments only. Clearly, in full 
analogy to (5.35), we have with <\\si{x,y) large enough 

C-^gJ [x^y]< df^ [x^y]< Cdf^ [x ^ y] (5.38) 
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for some C — C(/3, 5) > 1 independent of x, y. 
5.2 Random walk representation 

The quantity ^ y] admits a particularly useful interpretation in terms of a killed 

continuum random walk in environment with random obstacles. To see it consider a 
continuous-time random walk Zt-M2{x,s) '■= Zt in independent of ^^^1 and governed by 
the following dynamics 

• between critical events specified below move in a constant direction with speed 1, 

• with intensity given by 4 times the covered length element update the movement 
direction, choosing the angle (j) € (0, 27r) between the old and new direction according 
to the density | sin(0)|/4. 

The starting point and the initial velocity direction for Zf are chosen by taking a straight 
line I crossing B2(a;, S) according to the measure | / nB2(a;, 5) ^ 0}). The starting 

point of Zx is now taken to be one of the intersection points of I with (5B2(a;, 5), each picked 
with probability 1/2, while the initial velocity vector lies on I pointing outwards Miix^S). 
Let Zt = Zt-n2{x,5) be the random walk Zt killed whenever hitting its past trajectory. 
A crucial observation is that the probability element of the walk Zt containing a given 
polygonal path 9 e (^^^y;^ as its initial subpath is exactly 

k 

2M0 I lnl(.,S) ^ 0}) «^P(-41engthW)ndMiN), (5.39) 

where ei, . . . , are the segments of 9 while /[cj] stands for the straight line determined by 
Cj. Indeed, the prefactor [2fi{{l \ I fl M2{x,6) ^ 0})]""^ comes from the choice of the initial 
segment of Zt [the line on which it lies and one out of two equiprobable directions, whence 
the extra 2~^] while for the remaining segments we use the fact that, for any given straight 
line lo, we have ^{{l \ I (1 Iq E di, Z(/,/o) G d(f)}) = \ sin (f)\did(f) with di standing for the 
length element on Iq and with Z(/o, denoting the angle between / and Iq, see Proposition 
3.1 in [2] as well as the argument justifying the dynamic representation of the Arak in 
Section 4 ibidem and the proof of Lemma 1 in Schreiber [15] . Note that the direction update 
intensity for Zt was set to 4 to cancel out with the normalising constant J^^ \ sin0|(i0 = 4 
in the density | sin0|/4 for the new angle choice. Clearly, the formula (5.39) is also valid 
for Zt replaced by Zt since the paths in C^^^'^ are by definition self-avoiding. Thus, taking 
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into account that, by standard integral geometry, //({Z | I fl B2(x, S) ^ 0}) = 27r(5 and 
recalling (1.11) and (1.12) we rewrite (5.39) as ^[©^*^^''']Pl(d^). Consequently, recalling 
the definition of C^'^'^'^^ and using (5.39) we come to 

Lemma 5.1 For each C C C^^2/;<5 following representation formula is valid for the 
value o/[e^^2'''5][2](C) 

[e-^y''f\C) = 47r5Ecard{t > | Z,- e (Zi),>o n,^ dM^iy, 6), Z[o,i] e C}, 

where (Zt)j>o flin 9B2(|/, 5) stands for the collection of entry points of into 182(1/, 5), with 
exit points not taken into account. 

A simple yet useful conclusion of Lemma 5.1 is that, denoting by Z^^ — Z^^^^^^ the random 
walk Zt killed at rate /3 and, in addition, killed whenever hitting its past trajectory, we 
have for /3 > 2 

[e^-y''f\C) = 47r5Ecard{t > | zf''^ e {zf-%>o n,^ dM^iy, S), z]^-'^ e C}. 

(5.40) 

Consequently, writing now Z^^ — Z^^^^^,^ for the random walk Z^ killed at rate (3 and, in 
addition, killed whenever hitting its past trajectory or a contour of in view of Lemma 
5.1 and (5.40) the definition (1.13) yields 

T^^j[x ^ y] = 47r5Ecard[(zf -'l)t>o dB^iy, 5)]. (5.41) 

A similar representation can be provided for Tj^^^^y by additionally killing the random 
walk whenever it hits dD. A corresponding representation for Tj-^^ and T^^^j^^ can also be 
given, yet we omit it because it is unessential for our further purposes and involves certain 
technicalities due the fixed last segment. 

5.3 Finite volume approximations 

The following lemma shows that T^g^lx <-> y] is well approximated by Tj^g^jy^[x y] for 
sufficiently large domains D. We write Il{x <-> y; S) for the square of sidelength 2S + 
dist(x, y) with one pair of its sides parallel and equidistant to [x, y] and with the remaining 
two sides at the distance 5 from x and y respectively, perpendicular to [x,y]. 
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Lemma 5.2 For each sufficiently large (5 > 2 there exists a constant C — C{(3,5) > 
such that 

whenever dist{x,y) is large enough. 

Proof The relation Tj^^jj^^^^.gy^lx ^ < Tj^s)\-^ ^ ^] obvious and only the remaining 
inequality Tj^g^ [x y] < CTj^^^jj^^^y.^^^ [x ^ y] requires verification. In view of the random 
walk representation (5.41) it will follow as soon as we show that 

Pm(^^y;S)) - 1/] > C-'P^^j [x ^ y] (5.42) 

for some C > 0, where 

^ y] P {{Zf~%>0 reaches dM^iy, S)^ 

and 




y] := P ({zf ^')t>o reaches dM2{y,6) before hitting dD 



Indeed, it is easily argued that upon hitting dE2{y, S) once, the random walk Zj. is un- 
likely to hit it too many more times and, consequently, the expectation on the right-hand 
side of (5.41) is bounded above and below by some constant multiplicities of the probability 
on the right-hand side of (5.42), the same observation holds for the corresponding repre- 
sentation of the finite- volume quantity Tj^s^jjf^x^y5))\p^ ^ ^l- omit the tedious technical 
details of this conceptually simple argument. 

The proof of (5.42) splits into two parts. First, denoting by Ri{x <-> y;S) the infinite 
strip between the lines determined by the sides of Il{x <-> y; 5) perpendicular to [x, y], we 
show that 

P(g [x^y]< CP([\ [x ^ y] (5.43) 

for some C > 0. Below it will be convenient to use the name x-line (resp. y-line) for the 
boundary line (side) of Ri{x ^ y;6) at the distance S from x (resp. y), perpendicular 
to [x,y]. Next, writing R2{x y;S) for the infinite strip contained between the lines 
determined by the sides of Il{x -(-^ y; S) parallel to [x, y], we show that 
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as dist{x,y) — > oo. Write 

Combining this with (5.43) and (5.44) yields (5.42) as required for completing the proof of 
the lemma. 

To proceed with the verification of (5.43), on the event that the random walk Z^~'^^ 
reaches dM2{y, S) before being killed we decompose its trajectory into three subpaths 

• Cx^y;S '■= {zf '^^)[T^,Ty], where Ty is the first time 2jf~'^ hits the y-line while is the 
last time T^f hits the x-line before Ty. 

. C. (^f-'^)[o,r.], 

• '■= i^f '^^)t>Ty, 

with the additional convention that Tx ■— if does not reach the x-line and Ty :— +oo 

if zf~^^ does not reach the y-hne (we set respectively Q and (y := $ in these cases). 
On the a;-line we construct a double sequence {xi)i^^ of points with Xj+fc lying at the 
distance \k\5 from Xi (say above for A; > 0, below for k < 0) and with xq coinciding with 
the intersection point of the x-line and the line extending [x,y]. The sequence {yi)i(zi on 
the y-line is constructed in the same way and ordered in the same direction as {xi). Let x 
denote the point in {xi)i(zi which lies the closest to zj-f if r^; > and x := x otherwise. 
Likewise, let y be the point in {yi)iez lying the closest to zl-^ ^' if Ty < +oo and y '■= y 
otherwise. With this notation it is easily seen that 

xe{x}Li{xi, i&} y&{y}ii{yj,3&} (5.45) 

where [x ^ x;y ^ y\ stands for the supremum over the possible realisations of C,x^y 
connecting ]B2(x, 5) with B2(y, 5) of the conditional probability, given C,x^y, that the random 
walk if "^^ connects B2(x, 5) to B2(x, 5) and B2(y, 5) to B2(y, 5). Since is killed with 
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the constant rate (5 — 2 > 0, for arbitrarily small e we can find A = A(e) such that, uniformly 
over X, y with dist(x, y) large enough, 

Xi, dist{xi,x)>X Vj Xi yj, dist(yj,y)>\ (5.46) 

Putting (5.45) and (5.46) together yields 

Xi,dist{xi,x}<\;yj,dist{yj,y)<\ (5-47) 

For dist{x,y) large enough the double sum in (5.47) can be bounded above by some con- 
stant C[\] times Pls-B.i{x^y5))\-^ ^ v] because each path of connecting dM2{xi,d) to 
dM2{yj,d) in Ri{x ^ y;6) with dist(x, Xj) < A and dist{y,yj) < A can be modified into a 
path connecting dM2{x, S) to dM2{y, S) in Ri{x <-> y; S) by an appropriate surgery between 
X and Xi and between y and yj at a probability cost depending only on A. It seems natural 
to expect that the supremum sup^,. j^^. P(s^iij^{x^yS))\-^i ^ ^-dniits an upper bound very 
close to Pl;J [ X > y], because dist(xj, y^) > dist(x, y) for all x^, yj. While we are not able to 
establish such a bound, we easily show that there exists a positive constant C with 

P(ski^^y;S)) ^ < ^ 2/] (5-48) 

uniformly in x,y with dist(a;, y) large enough. Indeed, this is done much along the same 
lines as in the considerations leading to (5.45) and (5.47), so we only sketch the argument 
omitting technical details. We split each path of the random walk Z^^^^^l.. s) connecting 
dM2{xi,S) to dM2{yj,S) into two subpaths: the initial subpath (i connecting dM2{xi,6) to 
some dM2{z,6), z G (5Z^ with | dist(a;j, 2;) — dist(a;,y)| < 6 (in fact, z can be chosen as 
the 5Z^-lattice point closest to the point where the random walk -^j!^"^],, 5) first gets at the 
distance dist(a;, y) away from Xi) and the remaining subpath (2- Integrating over (i for 
fixed z yields a value bounded above by a constant multiplicity of Pjf^[x ^ y] with this 
prefactor (arbitrarily close to 1 for 6 small enough) due to the fact that dist(a;j, z) differs 
slightly from dist(a;,y). Integrating over (2 conditioned on ^1 and summing over z yields 
only a constant prefactor - the sum of integrals converges due to the constant killing rate 
/3 — 2 > along (2- This proves (5.48). Combining now (5.47) with (5.48) and with the 
discussion directly following (5.47) we obtain 

pI^][x ^y]< C[A]i^[\(_^^,))[a: ^ y] + eC'P^^'^x - y]. (5.49) 
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Choosing e small enough so that eC < 1 (recall that C does not depend on A) completes 
the proof of (5.43). 

To establish (5.44) we denote by v^y the unit vector pointing from x to y, i.e. Vxy '■— 
{y — x)/ dist(x, I/), and for small > we consider the event <-> y\ that 



9182(3;, 5) to 9182(1/, 5) cannot exceed [dist(a;, y) + 25]/[l — rj\. Since ^[^^ is stochastically 
dominated by the Poisson contour process Vqw, as follows by the graphical construction 
of Section 1.2, we conclude that there exists k > such that for all j3 large enough 



Indeed, to see it we: 

• Split the strip Ri{x <-> y]5) with equidistant straight lines perpendicular to [x,y\ 
into 6(dist(a;, I/)) equal-sized strips. 

• Construct a path of the random walk Z connecting dB>2{x,5) to dE>2{ii,5) and such 
that the scalar product of v^y and the current velocity vector of Zt falls into [1 — 77, 1] 
for all time moments before 9B2(y,(^) is reached. This is done by constructing and 
patching together subpaths of Z crossing individual strips, at a constant probability 
cost per strip. 

• Use the stochastic domination of J\}^^^ by Vqw to conclude that the probability that 
the so constructed path of Zt avoids ^['^l is bounded below by exp(— 0(dist(a;, 7/))). 

• Check for survival of the so constructed path under /3 — 2-kiUing, which yields a 
probability prefactor bounded below by exp(— ^5^[dist(x, y) + 25\). 

By the definition of ^'ff.^y for d\si{x,y) large enough this procedure allows us to bound 
below the probability of this event by exp (fr^ + K)[dist(x, y) + 25] j for some /t > 0, 
as required. Since Vqip] stochastically dominates "Pqi/s'i for {3' > P, this technique works 
uniformly in j3 large enough. We omit tedious technical details of this standard argument. 
To proceed, define the event TZ\f^ [x ^ y] that 



• The random walk {zf ^')t>o reaches dM2{y,S), 

• The scalar product of v^y and the current velocity vector of zf~^^ is in [1 — 77, 1] for 
all time moments t >0 before dM2{y, S) is reached. 

Observe that on the event Sl^K[x ^ y] the total length of the path of connecting 




(5.50) 
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• The random walk (Z| ^)t>o reaches dE2{y,S), 

• The random walk {zf ^')t>o hits dR2{x ^ y]S) before reaching dM2{y,S), 

and observe that on T^-'^jfx ^ y] the length of the path connecting d'E2{x, 5) to dM2{y, 5) 
has to exceed •\/5/4dist(x, y) — 25 and, hence, 

P (n^l^ix ^ 7/]) < exp - 2)[y574dist(a;,7/) - 25]) . (5.51) 

Noting that P (^S^^^^^-^[x ^ y]) < ^ ^] sufficiently small rj and putting 

(5.50) together with (5.51) we see that, for (3 large enough. 

This yields (5.44) and hence the required relation (5.42). The proof of Lemma 5.2 is 
complete. □ 

5.4 Existence and finiteness of surface tension 

In this subsection we use the preceding results of this section to show that 

Lemma 5.3 The limit defining the surface tension functional in (1-15) exists, is finite 
and strictly positive. 

Proof The main work has already been done in Lemma 5.2. In view of the relation 
s^PA>25'^i^' < oo as easily deduced from (5.50), the required existence of the limit in (1.15) 
will follow by a standard almost-subadditivity argument once we establish the following 
auxiliary lemma. 

Lemma 5.4 For £'[Ai, A2] := (log[Ai + A2])^ and Ai, A2 large enough we have 
(Ai + A2 + D[K A2])r^U+i.[A.A.] < AirJJl + X.rf} + 0{D[\,, A^])). 

Proof of Lemma 5.4 For fixed 5 > consider disjoint squares IIi := 11(0 <-> Aie^;; 5) 
and 112 '■— n((Ai + D[Ai, A2])e-r <-> (Ai + A2 + -D[Ai, A2])ea;; 5), separated by a moat of width 
-D[Ai, A2], and observe that, since -D[Ai, A2] ^ log[Ai + A2], by the decoupling Lemma 4.1 
applied to Ui :— Hi fl 9, U2 :— 112 H 9, with 9 standing for the integrand polygonal path in 
the definition (1.13) of rf^'[- ^ •], it follows that 

41[0^(Ai + A2 + i?[Ai,A2])e,] > 
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g-o(D[Ai,A.]) ^M^^^ ^ Aie,]7;^\)[(Ai + D[X^, AsDe, ^ (Ai + A2 + D[\^, A2])e,], 

(5.52) 

with the prefactor e~'^i^l^i'^^]) due to patching together pairs of paths 9i in C^^^^^""'^ and 
6*2 in C^^i+Dl>^u>^2])e,'^{Xi+X2+D[XuX2])e.;S^ ^^^^^ disjoint with ^l/^l, into paths 9 falhng into 
^o^(Ai+A2+-D[Ai,A2])ea:;(5 disjoint with A^^\ by constructing a path connecting and 62 across 
the moat of width D[Ai,A2] separating Hi and 112, according to a procedure completely 
analogous to that used in the argument leading to (5.50). Note that the fact that the 
patching procedure involves here conditioning on A^'^^ being disjoint with 9i and 62 does 
not affect this argument because the conditional graphical construction of the process 
A^l.Q^yjg^ guarantees that it is stochastically bounded by Vq^.q^^q^ and hence by Vqip] as 
used in the proof of (5.50). To proceed, apply Lemma 5.2 to conclude that the quantities 
^(5) [0 ^ ^i^x] and T^j[{Xi + D[Xi, AaDe^ ^ (Ai + A2 + L>[Ai, AaDe^] are bounded above by 
constant multiplicities of their respective finite volume counterparts T^}ji^-)[0 Aie^:] and 
T^}-^^^[(Xi + D[Xi,X2\)ex (Ai + A2 + -D[Ai, A2])ej;]. Combining this conclusion with (5.52) 
shows that T^^[0 ^ (Ai + A2 + L'[Ai, A2])e,] > exp(-0(L>[Ai, A2]))7;^^)^[0 ^ Xie,]T^j[{Xi + 
D[Xi, A2])ea; <-> (Ai + A2 + -D[Ai, A2])e^] for some C > 0, which completes the proof in view 
of the definition (1.14) of r^f , i = 1, 2. □ 

Completing the proof of Lemma 5.3 With the existence of the limit in (1.15) estab- 
lished we now easily conclude its strict positivity from the positivity of killing rate in the 
random walk representation (5.41) while the finiteness of r'^l follows by the probabihty 
lower bound (5.50). □ 

6 Skeleton estimates 

The purpose of this section is to provide coarse-graining estimates based on skeleton cal- 
culus. For a, 5 > 0, always assumed to satisfy L ^ a ^ S and to tend to 00 as L — > 00, 
by an (ct, 5)-skeleton in M2{L) we shall understand a collection (7i, Ei, I2, E2, . . . , Im, E^) 
of pairwise different points {skeleton vertices) in M2{L) fl Z^, with 7i, /2, ■ ■ ■ referred to as 
the initial points, Ei,E2, . . . as the corresponding endpoints and [/i, Ei], [I2, E2], ... as the 
skeleton segments, where the following is satisfied for alH = 1, . . . , m 

(SI) a-V2< dist(/i, Ei) <a + V2. 
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We say that a collection 7 of a-large polygonal contours is compatible with an {a, S)- 
skeleton E = (7i, Ei, . . . , 7^, E^n), write 7 ~ E, if the following holds for all i = 1, . . . , m 

(52) There exists a contour 6i ^ and points I^jEj G 6i such that dist(/j,/7) < 
^,dist(E,,E7) < ^ and dist(/i,a;) < a + ^ for aU a; e 9^11^, Ej] with 9i[I^ , E^] 
standing for the polygonal path from to E^ along 9^ (note that we do not require 
that 9i 7^ 9j for i 7^ j), 

(53) Either we have dist(/j, {/i, . . . , /j-i}) < a + 6 + \/2 or i is the smallest index with 
77 e for some 9i E ■j, 

(54) For each x e 7 we have dist(a;, {/i, . . . , 7^}) < 2q; + 5 + -\/2, 

(55) The polygonal paths 9i[I^ , Ej] are in a distance at least 5 away from each other. 

Roughly speaking, the motivation underlying this definition is the following. For two 
distant points x, y with dist(x, y) = fl{a) connected by a polygonal subpath of a contour we 
want to find a collection of approximately equal-sized segments [/j, Ei] of length a(l + o(l)), 
lying on this path and such that their overall length is at least dist(,x, y)(l+o(l)). Being only 
concerned with this total length condition, as ensured by (S3) stating that the distance 
between initial points is close to the single segment length, we do not require that these 
segments form themselves a connected polygonal path or that their ordering agree with 
the orientation of the path. On the other hand, we do impose an explicit lower bound 
(S5) for distance between polygonal subpaths crossing different segments, thus ensuring 
the applicability of the decoupling Lemma 4.1 in our further argument. It should be 
emphasised that this approach, considerably simplifying our argument in the sequel, can 
only work in an isometry invariant setting, as ours, where it is justified to look only at the 
total length of phase interfaces while ignoring their local directions, ordering etc. 

We say that a collection 7 of ct-large contours dominates an (a, 5) -skeleton E, write 
7 E, iff 7 contains a sub-family of contours 7' with 7' ~ E. Further, by the length of 
a skeleton E = (7i, Ei, . . . , 7^, E^n), denoted length(E), we understand the total length of 
skeleton segments Yl^i dist(7j, Ei). We write also A'"(E) for the total number of initial and 
endpoints in E. We say that a collection 7 of ct-large polygonal contours is well covered by 
an {a, 5)-skeleton E, write 7 oc E, if the following holds 

• 7 ~ E, 
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• E maximises length(E) among skeletons compatible with 7. 



For an {a, 5)-skcleton E we consider the corresponding black phase area, denoted in the 
sequel by Area(E), and given by 



In other words, Area(E) is the supremum value of possible black phase area which can be 
enclosed by a collection 7 of a-large contours well covered by E. We note that for some E 
there may be no such 7 in which case we put by convention Area(E) := 0. 

In the sequel when no ambiguity occurs we will often use the ~, oc notation for 
contour collection containing also a-small contours, in which case we shall always mean 
that the approriate relation holds for the corresponding sub-ensemble of a-large contours. 

Lemma 6.1 For a collection 7 of a-large contours there exists a compatible (a, S)-skeleton 



61 e 7, set 77 to be the point of 61 minimising the distance to Ii and let be the first point 
(say in clockwise order) on 9i at the distance a from Ij (note that the distance considered 
here and below is the usual Euclidean distance and not the distance along the contour ^1 !). 
Set El to be the point of B2(L)nZ'^ which hes the closest to Ej. The conditions (S1),(S2) 
for i = 1 is now easily verified. Further, if existing, choose I2 to be the point minimising 
the distance to 9i[Ii, Ej] (with ties broken in an arbitrary way) among the points I] in 7 
with the property that there exists E] e 62 with dist{l2, E]) = a, dist(/2,a;) < a for all 
X e 92[l2:E2] and dist(6'i[/7, -E^]) ^2[-^2 5 -^2]) ^ ^) where 62 is the contour of 7 containing 
II Note that if 9i = 62 then dist(/^, ^7]) = 5 and hence dist(/7,/7) < a + 5. In 
case such and £"7 exist, we define I2 and E2 as the best approximations in n 182 (L) 
of I2 and E2 respectively, getting the required relations (S1),(S2),(S3),(S5) for i = 2. 
On the other hand, if such a pair [I], E^) fails to exist, we conclude that no point of 7 lies 
further than a + 5 away from ^i[/7) ^i]^ otherwise we could find /J and E^ with desired 
properties. In this case by (S2) we have dist(a;, /7) < 2a + 5 for all a; in 7 which yields 
(S4), and (S3),(S5) are obvious, whence (/i,£'i) is already an (a, 5)-skeleton compatible 




E. 



Proof Choose an initial point /i e fl B2(-L) at a distance less that l/\/2 from some 



with 7. 
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We proceed inductively with this construction, adding new pairs (/j+i, -Bj+i) obtained 
as the best lattice approximations of {Ij^i, E'^^^) with /^fi arising as the point minimising 
the distance to [jj^-9j[IJ,Ej] among the points /^fi ^ ^i+i ^ 7 which there exists 
-E'7+i e 6i+i with dist(J2^i, -E^^+i) = a, dist(/j^i, x) < a for all x e 6'i+i[/2!^^, £"7+1] > such 
that dist(^i+i E7+J, U^.<. > 5. We note that if Oi+i = Oj for some 3 < i then 

dist(/7+i, [jj^iej[I], E]]) = 5 and hence dist(/7+i, {I^, . . . , Ij}) < a + 5. The construction 
terminates when no further pair can be found, and it is easily verified as in the argument 
above that the resulting collection (Ii, Ei, I2, E2, ...) is an (a, 5)-skeleton compatible with 
7. The proof is complete. □ 

Recalling that, by the definition, skeletons have their vertices pairwise different and 
belonging to B2(L) fl and hence their number is finite, we obtain the following corollary 
as an immediate conclusion of Lemma 6.1. 

Corollary 6.2 Each finite collection 7 of a-large contours can be well covered by some 
{a, S)-skeleton S. 

A particular feature of the notion of skeleton as introduced in this section is that if 
two polygonal subpaths of some contours go very close to each other, it may happen that 
only one of these subpaths will contribute to the total length of a well-covering skeleton 
because of the requirement that subpaths going along the segments of the skeleton keep 
distance at least 5 from each other as imposed in (S5) above. However, this does not lead 
to problems in our further argument, since we are mainly concerned with minimising the 
skeleton length given the enclosed area, where collections consisting of multiple contours 
are outperformed by singleton ones. This is made formal in the isoperimetric lemma below. 

Lemma 6.3 Assume that A <^ a^. Then for each {a, 5)- skeleton H in'E2{L) with AreaCE) — 
A, A e [0,7rL^], we have 

length(E) > 2V7r]4[l - 0{S/a)] - 0{a). 

Proof Below, we restrict our attention to skeletons S with length(E) < 27rL, since oth- 
erwise our assertion is obvious. 

Pick some collection of a-large contours 7* oc S with Area (black [(Jg»g^» 6*]) — Area(E) — 
0(1) = A — 0(1), and observe that, by the definition of the relation oc, to prove the lemma 
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it is enough to construct E* with 7* ~ E* and such that, for A:$> a'^, 

length(E*) > A^{A) (6.53) 

for a non- increasing function A 1— > ip{A) — ip{A;a,5) with A 1— > Aip{A) non-decreasing, 
satisfying 

iPiA) ^2y/^[l-0{S/a)]-0{a/A) (6.54) 

(note that the statement of the lemma triviahses for A — 0(q;^)). Without loss of gener- 
ality we can and do assume that 7* contains no nested contours, for otherwise we could 
simply remove the internal contours increasing the area enclosed by 7*, proceed with the 
construction below for the so reduced 7* obtaining S* of required length, and then con- 
struct some additional skeleton segments for the internal contours and add them to S* thus 
increasing its length even further. We also assume that 7* contains only contours for which 
dist(^^*, |J^*g^*^|g,| 5*) > 64a. This does not result in loss of generality because finding a 
sub-collection 7* of contours satisfying this condition and such that all other contours of 
7* are contained in 64Q;-ncighbourhood of IJ7*, and then constructing S* for 7*, we see 
that the total area enclosed by the contours in 7* \7* is of order 0(q; length(S*)), whence 
by (6.53) for S* we get length(S*) > [A-0(alength(S*))]^(A-0(alength(S*)) and con- 
sequently, by (6.54), length(E*) > lcngth(S*) > Aip^A) — 0{a) provided Q;length(E*) = 
0(^4). The remaining case length(E*) = fl{A/a) is easily handled directly, by considering 
subcases A — 0{a'^) and A » cc^. 

The proof of existence of E* satisfying (6.53) goes by induction with respect to the num- 
ber n* of contours in 7* = {O^, ... ,^**}, assumed to be ordered by decreasing enclosed 
area. For n* — 1 the assertion follows immediately by standard isoperimetric argument: 
note that the correcting term AO{a/A) — 0{a) coming to the RHS of (6.53) when substi- 
tuting (6.54) is due to the admissible distance ©(a) between a skeleton and a compatible 
polygonal path [see (S4)], while the prefactor l — 0{S/a) there comes from the fact that, in 
the single contour case, the distance between the initial point of a given skeleton segment 
and the set of preceding initial points may exceed the length of the segment by at most 
5 + 2y/2 [see (SI), (S3)], which is fraction 0{6/a) of the segment length. 

To proceed, take n* > 1. We split our argument into three possible cases. 

Case 1: Say that a point x e is a-seen from a contour 9* G 7* iff dist(a;, 6**) < Aa. 
Assume that the total length of the set seen(^**, a) of all such points exceeds 
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16a for some i < n* and recall that, as assumed above, there exists x E with 
dist(x, Uj<n* ^j) ^ Patching and 9* together with additional polygonal 

paths at two extreme points xi, X2 of seen(^**, 9*; a) and removing the internal parts 
of both contours between xi and X2, we replace 9^* and 9* by a single contour 9"^, 
which can be made disjoint with all remaining contours 9*, j i,j ^ n* . Denote by 
7^ the contour collection resulting from 7* by replacing 6'** and 9* by 9*^ and possibly 
removing some further contours which would become nested due to this replacement. 
It is easily seen that, by our assumptions above, any skeleton ~ 7^ can be 
modified into S* ~ 7* with length(S*) > length(E^). Thus, the assertion (6.53) for 
7* will follow if we are able to find such with length(S^) > Aipi^A). However, this 
is ensured by the inductive hypothesis in view of the relation Area(black[[Jg.g^* 9*\) > 
Area(black[Ug.g^. 9*]). 

Case 2: Next, suppose that length(seen(^**, a)) < 16a for all i < n* and that A^* » 
a^, where A^* stands for the area enclosed by ^**. Recall in addition that there exists 
X e with dist(x, Uj<n* ^j) > 64q;. We construct an (a, (5)-skeleton E* as follows. 
Put 7!. := 7* \ {^**} and observe that the area enclosed by 7!. is ^4 — An* — o(l), 
which is due to the fact that there is no contour nesting in 7* as assumed above. We 
let T,*_ be an (a, 5)-skeleton such that 

length(El) > - An*]ij{A - A^*) > [A - An*]^{A), (6.55) 

with its existence guaranteed by the inductive hypothesis [note that A — A„* ^ a^ 
since the contours are ordered by decreasing enclosed area] . The skeleton T,*_ can be 
extended to a skeleton E* compatible with 7* by the procedure described in the proof 
of Lemma 6.1. Denoting by S* :— E* \ El the collection of newly added segments we 
see by our assumptions for Case 2 that S* can be in its turn extended to an {a, S)- 
skeleton S* compatible with {6***} by adding at most 0{n*) new segments covering 
seen (6***, 6**; a). Thus, using isoperimetric argument, as applied for the case n* — 1 
above, we are led to 

length(5*) > A*V(^n«) - 0{n*). (6.56) 

Recall now that the contours 9^,92, ■■■ are ordered by decreasing enclosed area, 
whence 

An' < A/n*. (6.57) 
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Using (6.57) to rewrite (6.56) aslength(5*) > An*[ip{An*)-0{n*/An*)] > An*[tp{An*) 
—0{A/A'^t)] and then applying (6.54), noting that An* < A/2 by (6.57) and resorting 
to standard calculus in order to check that, for a large enough, we have tlj{An*)—0{A/ 
Al^) > il;{A) for An* > we conclude from (6.56) that length(5*) > An*i}{A) 
for An* ^ A^^^. On the other hand, the trivial bound length(5*) > a is easily seen to 
yield length (5*) > An*ip{A) whenever An* <^ ay/A. Since we assumed that A <^ 
in the statement of the lemma, we get a\/A 3> A"^/^ which leads to 

length(5*) > An*ip{A) (6.58) 

for all An* within range of (6.57). Combining (6.58) with (6.55) and recalling that 
length(E*) = length (El ) +length(«S*) yields the required relation (6.53) for Case 2. 

Case 3: Assume now that A^* = 0(«^). Then the required relation (6.53) can be obtained 
along the same lines as in Case 2 by recalling that dist(^**, IJ^-^^^, ^^*) > 64a and 
noting that putting such into 7* results in large added length to added area ratio, 
exceeding ip{A). The only reason for discussing this case separately is the technical 
fact that '0(^n*) is formally not defined for An* < Ca^ unless C is large enough. 

The proof is now complete by induction. □ 

With the concept of an (a, 5) -skeleton discussed above we are now in a position to 
proceed to the main result of this section. 

Lemma 6.4 With a — > cxd, 5 — >• cxo and L ^ a ^ 6 ^ ^ogL, we eventually have for each 
{a, S)-skeleton S in B2(I/) 

P {A^^ ^ E) < exp (-rf 1 length(E)) . 

Proof For a contour collection 7 >: E, 7 = {^1, . . . , 9k} we consider the partition E[7] = 
{»S[^i], . . . ,»S[^fe]} of E into disjoint sub-skeletons S[9j\ composed of segments [/, -E] with 
the corresponding points as given by (S2), lying on 9j. Note that some S[9j\ may 

be empty. Moreover, for a non-empty sub-skeleton 5 C E we write [S\ to denote the 
family of all contours 9 such that dist(/, 9) < l/y/2 and dist(£', 9) < l/V^ for all segments 
[I,E] e S. In particular, we always have 9j e [^[9^]] provided S[9j] ^ 0. In view of (1.6) 
or, equivalently, by the graphical construction of Section 1.2, we see that 

\j=l / j=l (6.59) 
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where the sum ranges over all possible partitions {Si, ... , 5^} of E and with the inequality 
rather than equality above due to the fact that we do not restrict the domain of integration 
to non-intersecting contours 9j and that more than one contour of A^^ might occur in [Sj], 
moreover it is not guaranteed that Sj — S[9j]. We fix a partition E = <Si U . . . U <Sfc and, 
to distinguish between vertices coming from different sub-skeletons Sj, j = 1, . . . ,k, we 
subscript skeleton vertices with the corresponding sub-skeleton names, writing Jj.^. and 
Ei-Sj- We also put 7 := {^^i, . . . ,0k}. Denote by Pi-j the polygonal subpath 9j[^\E^^] of 
the contour dj G [Sj] in the above integral, with /^^ and E^^ standing for the points of 
9j closest to li-s^ and Ei-^j respectively. Note that the points I,-^ and E^^ do not have to 
coincide with /^^^ and Ej^^ as specified by the correspondence (S2) implied by 7 ^ S; yet 
we clearly have dist{IJ.g,, 1^^) < a/2 and dist{Elg,, E^^) < \/2. The reason for introducing 
jf^ and E^^ rather than simply using 1^.^, and Ej^^ in their stead is to ensure measurable 
dependence of Pij on 9j. 

Observe that by condition (S5) the distance between different Pij does not fall below 
S. Given the collections {9) :— {9j) and (P) we consider the events 

i.,[9,] :=Mi^inP,, = 0}. 

Taking into account that Area(Pij ® B2(l)) = 0{a'^) (see (S2)) and using the decoupling 
Lemma 4.1 yields uniformly in (9) 

P ( n^*;j[^j] I = (1 + C'(exp(-C5)alength(S) logiV(E))). 

J i,j (6.60) 

To proceed, note that, by (6.59), 

{5i,...A}"^[5ilx-x[5*] Vm / i=l 

and hence, in view of (6.60), applying the rough bounds A^(E) = O(L^), a — 0{L) and 
length(E) = O(L^) we obtain 

P(^[^1^E)< Yl I \{^i^^■A0j])\{dQ^'\9,) 

{l + 0{exp{-CS)LHogL)). 
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For an endpoint Ei-Sj we write q{Ei.Sj) to denote the skeleton vertex liz-Sj or li'-Sj directly 
succeeding E^-Sj in clockwise order on 9j. Then, by the formulae (1.13) and (5.37) for 
T^^[- ^ •] and ^ •] respectively, in view of the definitions of T^^^l- •] and d^^^^[- •] 
as provided in Subsection 5.1, and by the definition (1.4) of the free contour measure and 
(1.11) of the free path measure we are led to 

. V ^v^^ ^v^^ / (6.61) 

(1 + 0{exp{-C6)LUogL))C^^^\ Ci > 0, 

where the inner sum ranges over all possible successor assignments q and where the extra 
factor < comes from integrating out the configuration of contours 6j within l/-\/2- 
neighbourhoods of li-Sj and Ei-s-, i — 1, . . . , which are subject to optimisation rather than 
integration in definitions of T^^^[- ^ •] and d^^[- <-> •]. Recall that dist(£'i;5^, ^■(E'i;^^)) > 
S — y/2 in view of (S5) and then use the random walk representation of Lemma 5.1 and 
(5.40) combined with (5.38) to conclude that Ja; <-> y] = exp(— $1(5)). Thus, taking 
into account that both the total number of possible partitions {Si, ... ,Sk} and the total 
number of possible successor assignments q are of order exp(0(A'"(E) log A'"(E))), in view 
of (1.14) the relation (6.61) combined with (5.35) gives us 

P(^[^] ^E) < J]exp(-T|^Mist(/,;5,,£^^;5j + 0(Ar(E)logAr(E)-Q(5Ar(E)))^ 

(l + 0(exp(-C5)LMogL)). 

Since, by definition, length(S) = '^■jdist{Ii.s.,Ei.s_^) and, moreover, exp(— C5)L'^ logL = 
o(l) and 5N{T,) ^ -/V(E) log A^(S) by the assumptions of the lemma, we conclude that 

P (^[^1 ^ E) < exp (-ri^l length(E)) 

for a,S,L large enough, as required. □ 

7 Lower bound 

Below, we provide a lower bound for the occurrence probabilities of large contours in A^'^^. 
This is complementary to the upper bounds obtained in the preceding Section 6. For 
a,S > and for a piecewise smooth closed curve a in we consider the event U[a;a] 
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that there exists a contour 9 e A^^ such that pij(cr, 9) < 2a with Ph{-: •) standing for the 
usual Hausdorff distance. The following lemma gives a lower bound for the probabihty of 
such event for a := = dM2{R). 

Lemma 7.1 With a — > oo, 5 — > oo, it! — > oo such that log R<^S<^a<^Rwe have 

P(W[§i(i?);a]) > exp {-27iRtI^^ - 0{6R/a)) . 

Proof Note that 

¥{U[Si{R);a]) > [ P n ^[^1 = 0) e[^l(d^) - P(W(>i)[§i(i?); a]), 

where W*^^^)[§i(i?); a] is the event that there exist at least two contours ^i, ^2, • ■ ■ in A^^^ 
such that p/f(§i(i?), 9i) < 2a, i = 1,2,... . Using the conditional graphical construction 
with forbidden regions we easily see that F(U^^^^[E>i{R); a]\U[Eii{R); a]) = o(l), whence 
(7.62) becomes 

P(W[§i(i?); a])> [ F{9n ^l^l = 0) e^f^\d9){l - o(l)). 

J {eec I (e,Si(i?))<2a} (753) 

To proceed, we partition the circle Si(-R) into disjoint segments Ei], i = 1, N{R; a, S) 
= Q{R/a) separated by spacings of length 5 and such that dist(/i,£'i) = a, i = 1, . . . , 
N{R;a,5). Denote by Ilj the square n(/j ^ Ei; 1/V2) as defined in the lines preceding 
Lemma 5.2. Clearly, Hi arc disjoint and dist(nj,nj) = 0(5) for i 7^ j. The integral 
in (7.63) can be bounded below by restricting the domain of integration to the family 
C[Ii, El, ...] of paths 6 such that, for all i = 1,... ,N{R;a,6), 9 contains a subpath 9i 
connecting (9B2(/j, l/\/2) to dM2{Ei, 1/V2) within Ilj. Using the decoupling Lemma 4.1, 
with Ui :— 9 n Hj there, we can factorize the integral 

F {9 n A^'^ = ^) e^'^\d9) 

{eeC[h,Ei,...] I ph{0Mo,R))<M 
into the product of T^ij^.^.^i^i E^], i = 1, . . . , N(R; a, 6) with a prefactor 
(l + 0(Rex.p(-Ci5)\ogN(R;a,5)))exp(0(5N(R;a,5))), Ci > 0, 

where (1 + 0{Rexp{— C 16) log N[R; a, S))) is the factorization correction from Lemma 
4.1 while exp{0{SN{R; a, S))) comes from patching the contour 9 by joining together the 
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subpaths 6i passing through adjacent (5-distant squares Hj so as to keep the resulting 
path within distance 2a from see the discussion of (5.52) and (5.50) above. Since 

Rexp{—CiS) log N{R; a, S) — o(l), we obtain 

N{R;a,S) 

F{U[S,{R)- a]) > exp{0{6N{R; a, 6))) J] ^(vV2;n,)[^^ ^ ^^]- 

i=l 

Applying Lemma 5.2 we conclude that 

N{R;a,S) 

P(W[§i(i?); a]) > OiC^^''-'''''^) exp{0{6N{R; a, S))) J] ^(VV2)[^^ ^ ^^]' > 

1=1 

Observing that N{R;a,d) = Q{R/a) completes the proof in view of the definition (1.14) 

8 Proof of the main theorem 

Throughout this proof we shall put 

a = a[L] := VZlogL and 5 = 5[L] := (logL)^. (8.64) 

As in the classical DKS theory, our argument below uses the decomposition of the contour 
ensemble A^^ fl B2(-L) into the collection ha-L ■— La;L [A^^] of a- large contours and the 
remaining family of ct-small contours, and it relies on an application of the skeleton bounds 
of Section 6 and complementary estimates of Section 7, combined with the use of moderate 
deviation results of Section 3. 

8.1 Lower bound for (1.17) 

In order to prove (1-17) we establish first the lower bound 



P (Ml(^[^]) > M{/3]7:L^ + aL\ N[a] L] holds) > exp | /^^Lt|^1 + 0{a] 



" 7(8.65) 



To show it, put 
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for some constant C large enough so that 

P {Ml{A^^) > M[(3]7rL^ + aL^ \ U[Si{R);a]) > 1/2, (8.67) 

with the event U[E>i{R); a], indicating the existence of a contour 9 of A^'^^ with p/f(^, 

< 2a, defined as in Section 7. Clearly, R < L for L large enough because a < 27r| M[/5]|. 

To see that the required choice of C in (8.66) is indeed possible note first that 

P (Mi(^[^l) < M[f3]nL^ + aL^ \ U[Si{R);a]) < 

Then use (3.20) to conclude that, for ph{9,Si{R)) < 2a, 

EMl (a^I^ ue^ = M[/5]7rL2 + aL^ + 47r| M[i3]\CRa + 0{La), 

which can be made larger than M[/3]7rL^ + aL^ by a term of order ©(La) under appropriate 
choice of C. In view of Corollary 3.3 this makes the integrand probability P(Mi(^^2.g U 
9) < M[/3]7rL2 + aL^) arbitrarily close to 0, uniformly in 9 with ph{^i{R),9) < 2a. 
In particular, (8.67) is seen to hold under such choice of C, as required. To proceed, 
observe that the conditional probability F{Af'^[a] L]\U[§>i{R);a]), being bounded above 
by [P(W[§i(;?); a])]-' /^,^,^ p„ieMR))<2a} ^(-^1"; ^] holds for A%)e^^KdO), tends to as 
L — s> oo by the results of Lemma 2.1 in Section 2 specialised for ^^ig- Thus, we conclude 
from (8.67) that for sufficiently large L 

P {Ml{A^^) > M[p]7rL^ + aL^, M[a; L] holds | U[Si{R); a]) > 1/4. 

The required relation (8.65) follows now by Lemma 7.1 in view of (8.66). 

8.2 Upper bound for (1.17) 

To complete the proof of (1-17) we shall establish the following upper bound, complemen- 
tary to (8.65), 

P (Mi(^[^l) > M\l3]nL' + aL', M[a; L] holds) < exp ( -J ^^Lr^^' + 0(a) ) . 

V V I^^^JI 7(8.68) 
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To this end use the exponential tightness bound in Lemma 2.1 to get for some Ci = Ci(a) 



P {Ml{A^^) > M[/3]7rL2 + aL^, X[a; L] holds) < 
P (Ml (^[^]) > M[/3]7rL2 + aL\ length(L,,.L) < dL, L] holds) + 0(exp(-C2L)). 



Applying Lemma 6.1 together with Corollary 6.2 we see that this probability is bounded 
above by 



^P (Mi(>ll^]) > M[p]nL^ + aL^, length(La;i) < CiL, L^^l oc E) + 0(exp(-C2L)), 



Note that we could restrict our attention to E C B2(L — [Aa — S — -\/2]) because of working 
on the event J\f[a,L], sec (S4). It should also be noted that any contour collection well 
covered by such E is completely contained in B2(L) for L and a[L] large enough. Under the 
imposed requirement that length(La;L) < CiL, the total length of skeleton segments in any 
E with F(ha-L oc E) > is also of order 0{L), whence the sum in (8.69) can be restricted 
only to skeletons of such length order. Observe now that the number of such skeletons 
is of order exp(0(Q;~^L log L)) which, by (8.64), is exp(0(\/L))- Indeed, constructing the 
skeleton segment after segment at each step we have at most 0{L^) possibilities of choosing 
new initial/end point. However, the total number of such steps, coinciding with twice the 
number of segments, is at most of order 0{L/a), because, as stated above, we only consider 
skeletons E with length(E) = 0{L) and the length of a single segment is close to a, see 
(SI). We put this statement as a remark for further reference 

Remark 8.1 The number of {a, 6) -skeletons mB2(L) wz^/i length(E) = 0{L) is of order 
exp{0{a~^LhgL)). 

Consequently, by (8.69), in order to estabhsh (8.68) it is enough to show that 




s 



(8.69) 



where the sum above is taken over all (a, 5)-skeletons E contained in B2 {L — [Aa — S — -\/2] ) . 



maxP {Ml{A^^) > M[(3]7rL^ + aL^, length(Lc,;L) < CiL, L^-l oc E) 




(8.70) 
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with the maximum taken over all {a, 5)-skeletons E satisfying the conditions specified above 
(i.e. contained in M2{L — [Aa — 6 — VO]) and with total length of order 0{L)). To proceed 
with the verification of (8.70) choose a skeleton Eq which achieves the above maximum. 
Putting 1/ :— Area(Eo) and A := length(Eo) we conclude by the isoperimetric Lemma 6.3 
that 

A > 2^[1 - 0{6/a)] - 0{a). 
Using that u — 0{LF') and that 5L/a = a we obtain 



\>2^/^-0{a). (8.71) 

To proceed, recall that M[/3] e (—1, 0) and observe that on the event Lq;l oc Eq we get by 
(3.20) and (S4) 

EMi(^gp"^(^^ U 7) < M[/3](7rL2 - v) - vM[l3] + 0{La). 

Thus, noting that the field A^^"^ conditioned on Lq,;^ = 7 coincides in law with .4.j^i!"''^^'''^''U7 
and recalling that A^(S) = 0(A/a), we conclude from Lemma 6.4 and Theorem 3.1 applied 
conditionally on 'La;L that the probability maximised in (8.70) is bounded above by 



exp {—T^^X) exp 



a 



(8.72) 



where A := [M[/?](7rL2-z/) -z/ M[/?] +0(L«)] - [Mf/^JvrL^ + aL^] = 2z/| M[/?]| -aL2 + 0(La) 
is the difference between the expected and actual (required) magnetisation on the event 
La;L oc So, and with v and A related by (8.71). Recalling that A = 0{L\a = y/L\ogL, 
S = (logL)^ and applying the lower bound (8.65) we see that the maximum in (8.72) has 
to be reached with A = 0(L^/^logL) = 0{La) and, consequently. 



2\M[P]\ 

whence, by (8.71), 



+ 0{La), (8.73) 



with the equality rather than inequality in the last formula due to (8.65). By Lemma 6.4 
this yields the required relation (8.70) and hence completes the proof of (8.68). 
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8.3 Existence of a large contour 

In view of the lower bound (8.65), the argument leading to (8.72) with the optimal 
skeleton Sq replaced by a generic skeleton S shows that, conditionally on the event 
{Ml (^['^1) > M[/3]7rL'^ + aL^, Af[a,L] holds}, with probability tending to 1 we can have 
La,L oc S only for those {a, 5)-skeletons S which satisfy (8.73) and (8.74) with i> = Area(E) 
and A = length(S). By the definition of the relation oc and by the proof of the isoperi- 
metric Lemma 6.3 this means that with conditional probability tending to 1 on the event 
{Ml {A^^^) > M[f3]7rL'^ + aL'^, J\f[a,L] holds} there exists at least one contour ^large of 
length + 0(a) and enclosing area 2\m[i3]\ + 0{La). In fact, we claim that for K 

large enough, conditionally on {M^ (^1^1) > M[/3]7rL^ + aL^, J\f[a, L] holds}, with proba- 
bility arbitrarily close to 1 the contour ^large is the only Xa-large contour of A^^ in M2{L). 
Indeed, for each E as above, i.e. satisfying (8.73) and (8.74), we have 

P {La,L oc E, A^^ contains more than one Xa-large contour in M2{L), J\f[a, L] holds ) 
< / P (lengthL,,L > Kaj e[^l(ci^iarge), 

^large 

where the integral ranges over 6'iargc in M2{L) of length L^I |^t|t| +0{a) and enclosing area 



aL 



2| 



f 0{La). Now, Remark 8.1 and Lemma 6.4 imply that the total mass of 
such ^large's is of Order exp ^ — jSpJI -^'''"^^ -\- 0{a)^ . Moreover, by Lemma 2.1 applied to 
'^R2-6»i integrand probability is uniformly of order 0(exp(— Xo;)). We now conclude 

our claim for K large enough in view of the lower bound (8.65). 

8.4 Uniqueness of the large contour, excluding intermediate con- 
tours 

It follows by the previous Subsection 8.3 that, conditionally on the event {Ml (^''^1) > 
]7rL^ + aL^, M[a,L\ holds}, with overwhelming probability there exists one large 



contour 6'iargc of length L. , + 0(a), enclosing phase area 2\m8]\ + 0(La), and this is 



the only i^a-large contour of ^''^^ hitting ]B2(L), with K large enough. Below, we argue 
that for sufficiently large Ciargci with overwhelming conditional probability, ^largc is in fact 
the unique Ciarge log L-large contour of ^'^1 hitting M2{L). The first step in this direction 
is showing in Lemma 8.2, similar to Lemma 4.2.4 in [12], that the phase of i^a-large 
contours adjusts very tightly to the micro-canonical constraint Mi(^l^]) > [M[/?]7r -|- a]L^ 
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and, roughly speaking, 'not much work is left for small contours'. Next, in Lemma 8.3 we 
use this knowledge to deduce the uniqueness of the large contour ^large and to exclude the 
presence of any other Ciaxge log large contours with overwhelming probability under the 
micro-canonical constraint. 

To proceed with the first of the afore-mentioned steps, we claim first that 

Lemma 8.2 With K as specified above we have 

P ([M[/5]7r + a]L2 - E (M,.(^[^]) | Lxaji) > L^'^\ 
Mi(^[^l) > [M[/3]7r + a]L^A^[Q;,L] holds) = o(l). 

Proof We set 

Applying Lemma 6.1 and Corollary 6.2 we get 

P ([M[/3]7r + a\L' - E (Mi(^[^l) | W^^^) > L^'\ Ml{A^^) > [M[(3\'k + a]L') < 
P (Mz,(>l[^l) > [M[/3]7r + a]L\ [M[/?]7r + a]L^ - E (Mi(>l[^]) | W^.^) > L^/^, 

E 

with the sum ranging over all (i^a, 5) -skeletons E contained in M2{L) and with oc' 
used as an indexed version of oc to denote the well-covering relation of {Ka, (5)-contours 
by (Xq;, 5)-skeletons. Use the exponential tightness results of Lemma 2.1 in Section 2 to 
conclude that, with arbitrarily large Ci and with C2 large enough, this sum can be bounded 
above by 

^ P (^[M[/3]7r + a]L^ - E {Ml{A^^) \ l^Ka^L) > L'l\ l.Ka;L ^""^'^ 

E, length(S)e[Ly^^-p,C2L] 

+ J] P (^Mi(^[^l) > [M[(3]ti + a]L2, L^^.^ ^"^^^^ + exp(-CiL). 

E, length{E)<L^^^-p (8-75) 

We proceed by showing that all consecutive terms in (8.75), for brevity denoted below 
by P11P2 and P3 respectively, are negligibly small compared to the probability P4 := 
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F{Ml{A^^) > [M[P]n + a\L'^,Af[a, L] holds). To begin with the first term Pi, use Remark 
8.1 to conclude that the number of summands in this sum is of order exp{0{a~^L log L)). 
Moreover, applying Lemma 6.4, noting that conditionally on 'LKa;L — 7 the field A^^ 
coincides in law with A^l'.^"''^^^^^ U 7 and using Theorem 3.1 conditionally on 'LKa;L, we 
uniformly bound above each summand of Pi by 



exp 




I exp 



A 



^4/3' 



L2 Ka 



exp - 



27ra [s] 
M[p]\ 



0(p) exp (-cL2/3) . 



RecaUing the definition of a = q;[L] = -\/Llog-L, p — L'^/^'^ and using the lower bound 
(8.65) of Theorem 1.2 with a replaced there by Ka, we conclude that 

Pi = o(P4). (8.76) 

To show that 

P2 = o(P4) (8.77) 
observe that, by isoperimetric Lemma 6.3, length(S) < L,J^^^^—p implies that Area(S) < 
- fi(Lp), whence, by (3.20), E {Ml{A^I^^ 1^k<.;l) < [M[l3]n + a]L^ - Vt{Lp) almost 



2|M[/3]| 



{Ka,&) 



surely on the event {Wa;L oc' S}. Consequently, by Theorem 3.1 applied conditionally 



on t^Ka;L, ©ach summand in P2 is bounded above by exp y—c ^ -£ j = exp (— c-^) . 
Using the lower bound (8.65) of Theorem 1.2 (with a replaced there by Ka) and recalling 
that ^ » L we obtain (8.77). Observing that, by the same lower bound (8.65), P3 — o(P4) 
provided Ci is chosen large enough, we complete the proof of the Lemma by combining 
(8.76) and (8.77). □ 



As announced above, our next statement will allow us to exclude with overwhelming 
conditional probability under the micro-canonical constraint the presence of Ciargelogi^- 
large contours different than ^large- 

Lemma 8.3 There exists a constant Ciarge > such that uniformly in collections 7 of 
Ka-large contours in B2(L) with length(7) < LlogL and in A < Lr'l^ we have 



P 



(^3;:^"'"'^'^^ contains a darge log L- large contour I (^^1;;^"'"^^'^^ U 7) > 

EM^ (^g::^"'"^^''^ u 7) + a) = 0(1). 
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Proof For brevity write 



..Ka ._ 
fJ'L,j - 



and let ^^[Ciarge, be the event that A^l'.^'^^^^^ contains no Ciargc log -^-large contours 
hitting M2{L) and S^lC^r^gc, L] the event that A^^^j^y"'^^^^^ contains no Ciargc log -^-large 
contours hitting 182 (L). From Corollary 3.5 it follows in particular that for each rj G 
[A, L*^/^ log L] there exists a unique value of the external magnetic field h[r], L] :— h[r], L, 7] = 
e(77/L2) such that 

E/.f;''^!"'^! = + (8.78) 

and, moreover, h[r],L\ increases with 77 given L, whence h[r),L\ e [/i~[L], /i"'"[L]] with 

h-[L] h[A,L]^e{A/L^) 

and 

h+[L] := h[L'^/'^ log L,L] = Q{L~^/HogL). 

For each L > we split the interval [A, L^/^ j^g^] ^^^^ Q{h+[L]L^/^ logL) = &{L'^/^ log^(L)) 
equal-sized subintervals [A = 770, ?7i), [r^i, 7^2), ■ ■ • of length Q{l/h'^[L]) — 6(L^/^/logL) 
each and we put hk^i '■— h ^'Uhl^+l^ x^j . Por each of the subintervals [rjk, rjk+i) write 

exp(-/lfc,i[E/xf;^ + rjk])r (^'^'^-^aarge, L] holds, //f^^''^'^'" e E//f_^ + 77fe+i)) 



exp(-/i,,4E/.g!^ + %+i])P e E/.g:^ + 77fc+i)) 

exp(/^,,^[77,+i - 77,])P (^£:'^'=^Maa.ge, - E/.f;''^^'-| < !?^±1-ZZ^^ 



< 



(8.79) 



At this point we claim that 
P 



I^I^Ka,^,,, _ i^i^^-^L^ < Vk+i_Vk}j ^ Q(i,-i/3iog-i L) (8.80) 

uniformly in L, 7, A, /c. Note that this is in fact a rather weak statement in the spirit of 
local central limit theorem (LCLT) and an LCLT could in principle be established for 
the polygonal Markov fields in its full strength much along the same lines as Lemma 
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2.4.1 in [12], with standard modifications due to the non-lattice nature of our setting. 
However, since we only need the weaker relation (8.80), we provide a much shorter argument 
speciahsed for this case. To this end, we subdivide the disk M2{L) into Q{L) equal-sized 
squares Qi,l,Q2,l, ... of side length ©(-s/Z), separated by moats of width log^(L). Now, 
in view of (1.10), the family of identically distributed random variables 



can be coupled with a sequence of i.i.d. copies Xj ^ of Xj ^ in the way that P(3jXj^£, 7^ 
Xi^L) = 0{L^exp{-c\og^{L))), c > 0. Indeed, 0{L'^ exp{-clog^{L))) is the order of the 
probability that the ancestor clans arising for different Qi^L in the graphical construction 
of Subsubsection 1.2.2 are not all pairwise disjoint. Write 



Yr 



[/3,/ife,J;Js:a,B2(L) 



U7) 



and note that 



(8.81) 



Further, observe that, in complete analogy with Theorem 3.1, 
F{\Yl - EYl\ > v^log'(L)) < exp ( -c 



[VL\og%L)f ^ VLlog^L) 



L\og\L) 



Ka 



— e 



-clog2(L) 



(8.82) 



Using the coupling of Xi^^ and Xj.i as discussed above, taking into account (8.81) and 
(8.82) and recalling that i^k+i — Vk = 0(i^^^^/ logL) ^ ^/L\og'^(L) we can now deduce 
the required relation (8.80) by the classical local central limit theorem applied for 
use e.g. Theorem 1 in Wey [17] with Al := L'^^^/\og{L) and Ml := L2/3/log^(L) there, 
with the assumption (HI) [central limit theorem for /U^"'^'"'^] there following by Theorem 
2.10.5 of [10] or Theorem 2.4.R5 and Section 5.3 in [9] with obvious modifications due to 
the continuum rather than lattice nature of our setting, and with the assumption (H2) in 
[17] satisfied in view of (3.25), (3.27) and by the relation Var[/x^"''**'^] = 6(L^) which can 
be established along the same lines as (3.29). 

Consequently, since hk,L['>lk+i — Vk] = 0{1) in (8.79), combining the relation (8.79) with 
(8.80) and taking into account that P(£^'*'='^[Ciarge, -^]) = exp(— n(CiargelogL)) uniformly 
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in A, L, /c, 7 for Ciarge large enough in view of Lemma 2.1 in Section 2, we conclude that, 
uniformly in A, L, /c, 7 

P(£:[a,,ge, e E//f^^ + [rjk, Vk+i)) = 0(1). (8.83) 

In view of (8.83) the assertion of our lemma will follow as soon as we show that 

p(/xf,:^ > E/xf;:^ + l'/' log li/x^ > E/xf;^; + a) = o(i) (8.84) 

uniformly in A, L, 7. To this end, use Theorem 3.1 to conclude that 

P(aiS > %S + log^) < exp(-cL2/3 log'(L)). (8.85) 

Next, apply Lemma 3.6 to get 

P(//f;^ > E/if;^ + A) > exp(-0([A + LlogLjVL^)) > exp(-0(L2/3)). 

which yields the required relation (8.84) when combined with (8.85). The proof of the 
lemma is hence complete. □ 

Recalling that conditionally on I^Ka;L — 7 the field J\^^ coincides in distribution with 
«4^i;;^°'''^^^^^ U7 and that, by the discussion in Subsection 8.3, 'hKa,L — {^large} with over- 
whelming probability under the micro-canonical constraint, and then combining Lemma 
8.2 with Lemma 8.3 applied conditionally on 7 = ^Ka,L-i shows that, conditionally on 
the event {M^ (^1^1) > M[/3]7rL^ oL^, H[a, L] holds}, with overwhelming probability 
^large IS the Only Cmax log iv-large contour of A^^^ hitting M2{L). This completes the present 
subsection of the proof. 

8.5 Localising the large contour 

It remains to show that the large contour 6'iarge satisfies 

But this follows immediately by specialising to our setting for ^^large the inequality (2.4.1) 
in Section 2.4 of Dobrushin, Kotecky & Shlosman [7] and combining it with (8.73) and 
(8.74). This completes the proof of Theorem 1.2. □ 
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9 Appendix 



Below, we discuss the dynamic representation and some further properties of the basic 
Arak process, see Arak & Surgaihs [2], Section 4 for the dynamic representation. For a 
fixed bounded open convex domain D we shall construct the basic Arak process A}) with 
free boundary conditions (unlike in (1.1) where empty boundary conditions are imposed). 

9.1 Dynamic construction of the basic Arak process 

We interpret the domain D as a set of time-space points {t, y) G D, with t referred to as the 
time coordinate and with y standing for the spatial coordinate of a particle at the time t. In 

this language, a straight line segment in D stands for a piece of the time-space trajectory 
of a freely moving particle. For a straight line / non-parallel to the spatial axis and crossing 
the domain D we define in the obvious way its entry point to D, in(/, D) e dD and its 
exit point out(/,Z}) G dD. 

We choose the time-space birth coordinates for the new particles according to a homo- 
geneous Poisson point process of intensity ti in D (interior birth sites) superposed with a 
Poisson point process on the boundary (boundary birth sites) with the intensity measure 

k{B) = Ecard{Z e A, in(/, D) e B}, B C dD. (9.86) 

Each interior birth site emits two particles, moving with initial velocities v' and v" chosen 
according to the joint distribution 

e{dv\dv") ■.= Ti-'\v' -v'\{l + v'\^'\l + v"^)-'"^dv'dv". (9.87) 

This can be shown to be equivalent to choosing the directions of the straight lines repre- 
senting the space-time trajectories of the emitted particles according to the distribution 
of the typical angle between two lines of A, see Sections 3 and 4 in [2] and the references 
therein. It is also easily seen that the value of angle G (0, tt) between these lines is 
distributed according to the density sin(0)/2. Each boundary birth site x G dD yields one 
particle with initial speed v determined according to the distribution dx{dv) identified by 
requiring that the direction of the line entering D at x and representing the time-space 
trajectory of the emitted particle be chosen according to the distribution of a straight line 
/ G A conditioned on the event {x = in(/, D)}. 

All the particles evolve independently in time according to the following rules. 
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(El) Between the critical moments listed below each particle moves freely with constant 
velocity so that dy — vdt, 

(E2) When a particle touches the boundary dD, it dies, 

(E3) In case of a collision of two particles (equal spatial coordinates y at some moment 
t with (t, y) e D), both of them die, 

(E4) The time evolution of the velocity Vt of an individual particle is given by a pure-jump 
Markov process so that 

F{vt+dt & du \ Vt = v) = q{v, du)dt 

for the transition kernel 

q{v,du) := \u - v\{l + u'^y^^'^dudt. 

It is worth noting that, in full analogy with the discussion following (9.87), the (sharp) 
angle between the straight lines representing the space-time trajectories of the particle 
before and after the velocity update is distributed according to the typical angle between 
two hues of A. 

It has been proven (see Lemma 4.1 in [2]) that with the above construction of the 
interacting particle system, the time-space trajectories traced by the evolving particles 
coincide in distribution with the Arak process A*d defined as in (1.1) with the family Fd 
of admissible polygonal configurations extended to allowing also for partial contours 
chopped off by the boundary, which amounts to admitting not only internal vertices of 
degree 2, as in (P2), but also boundary vertices of degree 1. 

9.2 Properties of the basic Arak process 

As already mentioned in the introductory section, and as shown in Arak & Surgailis [2] , the 
basic Arak process A*j^ enjoys a number of striking properties. The two-dimensional germ 
Markov property, stating that the conditional distribution of the field inside a bounded 
region with piecewise smooth boundary given the outside configuration only depends on the 
trace of this configuration on the boundary (intersection points and intersection directions) 
is an immediate consequence of the Gibbsian definition. Next important property is the 
consistency: for bounded open and convex Di and D2 with Di C D2 the restriction of A}) 
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to Di coincides in distribution with A*jjj^, see Theorem 4.1 ibidem. This immediately allows 
us to define the infinite volume Arak process A, which inherits the isometry invariance of 
the finite volume Gibbsian definition and which is a thermodynamic limit for A^^^ . By the 
results of Schreiber [14], this corresponds to the unique infinite- volume bounded-density 
stationary evolution of the particle system discussed in Subsection 9.1 above. Interestingly, 
the intersection of the Arak process A with any fixed straight fine is a Poisson point process 
of intensity 2, see [2], which gives us direct access to two-point correlation functions of A 
under the colouring as in Subsection 1.1. Moreover, the partition function for the Arak 
process can be explicitly evaluated: it is known that 

E exp(-21ength(5)) = exp(7r Area(D)), 

5erj,(AB) 

see Theorem 4.1 in [2] [note that the prefactor 2 exp(length(9D)), present in the quoted 
theorem, is absent here because we take the law of A rather than the unnormalised measure 
II* as the reference measure and, moreover, we do not sum over two different admissible 
black/white colourings of each polygonal configuration]. It should be emphasised that 
these exact results are only available for A and not for A^^\ /3 > 0. 

Interestingly, there exists a much broader class of consistent polygonal Markov fields 
admitting analogous dynamic representations, possibly enhanced to allow for vertices of 
higher degrees (3 and 4), see ibidem. The question of characterising the class of all polyg- 
onal Markov fields admitting dynamic representations is far from being trivial and falls 
beyond the scope of this article. A conjectured description of this class has been provided 
in Arak, Clifford & Surgailis [4] , where a very nice alternative point- rather than line-based 
representation of polygonal fields is also discussed. 
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